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I.  INTRODUCTION 


In  the  failure  analysis  of  composite  laminates,  one  of  the  most  serious 
problems  has  been  the  propagation  of  interlaminar  cracking,  commonly  known 
as  delamination.  This  mode  of  failure  is  a  major  cause  for  the  deterioration 
of  laminate  structural  properties,  including  its  strength,  stiffness,  re¬ 
liability  and  durability.  Although  delamination  occurs  initially  as  a  local 
damage,  it  can  propagate  rapidly  to  cause  catastrophic  laminate  failure, 
especially  under  compressive  and/or  fatigue  loads. 

The  exact  formation  and  growth  mechanisms  of  interlaminar  cracking  in 
laminates  are  not  well  understood.  The  general  belief  is  that  a  certain 
distribution  of  small  interface  flaws*  exists  in  laminates  prior  to  loading. 
These  are  thought  to  be  caused  by  factors  such  as  poor  fabrication,  improper 
handling,  residual  curing  stresses,  external  impact  damage,  environmental 
degradation,  etc.  Under  a  certain  critical  loading  condition,  some  interface 
flaws  would  grow  and  coalesce  with  each  other,  forming  a  single  crack  of  macro¬ 
scopic  proportion.**  At  this  dimensional  level,  such  an  event  would  constitute 
the  onset,  or  the  initiation  of  the  macro-crack,  the  delamination. 

Experimental  evidence  indicates  that  upon  the  initiation  of  a  delamina¬ 
tion,  the  growth  process  may  remain  stationary  in  some  laminates;  while 


*  The  size  of  the  interface  flaws  is  of  the  order  of  the  fiber  diameter. 

For  graphite-epoxy  composite  systems,  for  instance,  this  dimension  is 
about  4  -  5  um, 

**  Observable  interfacial  cracks  are  in  the  realm  of  the  ply  thickness,  which 
is  two  orders  of  magnitude  larger  than  the  size  of  fiber  diameter  in  most 
commercially  available  graphite-epoxy  systems. 


NADC-79056-60 


in  others  the  process  may  be  unstable  and  lead  to  catastrophic  growth.  Clearly, 
the  factors  that  influence  the  nature  of  delamination  initiation  and  growth 
are  complicated.  The  exact  cracking  mechanisms  may  ultimately  involve  intri¬ 
cate  material,  geometrical  and  loading  interactions 

Early  investigators  focused  mainly  on  experimental  studies  of  free  edge 
induced  delamination  in  laboratory  test  specimens  [1-6].  The  emphasis  was  to 
determine  the  effect  of  laminate  stacking  sequence  on  tensile  strength  and 
fatigue  life.  It  was  found  that  the  laminates  of  particular  stacking  sequences 
in  which  free-edge  delamination  occurred,  always  had  a  reduced  strength. 

The  major  cause  of  free  edge  delamination  is  the  tensile  component  (normal 
to  lamina  interface)  of  the  interlaminar  stresses  that  concentrate  near  the 
laminate  edge  boundaries.  The  sign  and  magnitude  of  the  interlaminar  stresses 
vary  considerably  with  the  stacking  sequence  of  the  laminate. 

A  considerable  amount  of  effort  in  the  past  decade  has  been  expended  to 
calculate  analytically  the  free  edge  stresses  in  multilayered  laminates  [7-19]. 
While  these  edge  stress  solutions  tend  to  explain  qualitatively  the  observed 
edge  effects,  they  are  still  inadequate  to  provide  a  quantitative  description 
of  the  delamination  failure  process,  including  initiation  and  growth  behavior. 

Since  the  delamination  cracking  mechanism  is  essentially  a  fracture 
problem,  a  description  of  the  cracking  process  requires  a  fracture  mechanics 
analysis.  In  a  recent  series  of  papers  [20-23],  Wang,  Crossman,  et  al. 
reported  on  a  fracture  mechanics  study  of  the  free  edge  delamination  problem 
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in  graphite-epoxv  laminates  under  simple  tension.  The  analytical 

approach  employed  is  based  on  the  classical  fracture  mechanics  concept  of 
strain  energy  release  rate  [24].  A  numerical  (finite  element)  routine  incor¬ 
porating  the  technique  of  crack-closure  [25]  is  developed  to  simulate  the 
delamination  cracking  process.*  The  method  proved  to  be  viable  in  describing 
the  delamination  process,  correlating  with  an  experimental  verification  effort 
which  was  conducted  concurrently  by  Wang,  et  al.  [23]. 

The  main  results  obtained  in  the  Wang/Crossman  study  are  briefly  dicussed 
below: 

a)  One  of  the  basic  assumptions  in  the  Wang/Crossman  approach  is  that  pre¬ 
existing  interfacial  micro-flaws  propagate  and  coalesce  into  a  macroscopic 
crack  whenever  a  certain  condition  is  reached  due  to  the  applied  load.  In 
order  to  bring  the  microscopic  effect  to  the  macroscopic  level,  the  concept 
of  "equivalent  edge  flaw"  is  introduced.  Namely,  the  effect  of  the  preexist¬ 
ing  micro-flaw  distribution  near  the  free  edge  on  the  initiation  mechanisms 
of  delamination,  is  macroscopically  equivalent  to  that  of  a  single,  large 
crack  of  size  ac.  ac  is  thus  an  inherent  material  property  of  a  given  laminate 
system. 

b)  As  is  postulated  in  the  classical  fracture  mechanics  [24],  it  is  assumed 
that  strain  energy  is  released,  and  is  converted  into  surface  energy  whenever 
delamination  cracking  is  induced  by  the  applied  load.  The  rate  of  the  avail¬ 
able  energy  release  per  newly  created  cracking  surface  area,  G,  is  the  driving 
force  which  extends  the  crack  further.  Of  course,  the  capability  of  the 


*  The  numerical  technique  was  first  used  by  Rybicki  and  Kanninen  [26]  to 
calculate  the  stress  intensity  factor  by  the  crack-closure  procedure. 
Rybicki,  et  al.  [27]  applied  the  method  to  calculate  the  energy  release 
rate  in  a  stable  delamination  cracking  problem. 
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material  to  convert  the  strain  energy  into  surface  energy  represents  the 
resistance  (toughness)  against  further  crack  extension.  The  latter  is  the 
well  known  critical  energy  release  rate,  denoted  as  G£.  The  quantities 
and  a^  are  the  only  two  material  parameters  appearing  in  the  predictive 
delamination  growth  model  suggested  by  Wang  and  Crossman  [20]. 
c)  The  behavior  of  the  computed  energy  release  rate,  G,  depends  not  only 
on  the  applied  load,  but  also  on  the  interface  location  of  the  delamina¬ 
tion  and  other  geometrical  parameters.  In  particular,  the  amount  of  strain 
energy  available  for  release  is  controlled  largely  by  the  laminate  stacking 
sequence,  the  lamina  geometry  (thickness  of  lamina)  and  the  overall  laminate 
thickness  relative  to  the  size  of  the  delamination.  In  short,  G  is  strongly 
influenced  by  the  lamination  structural  dimensions  at  the  ply  level. 

To  summarize  the  above  development,  the  simple  case  of  free  edge  delami¬ 
nation  in  laminates  subjected  to  static  tension  has  been  experimentally  studied 
and  analytically  simulated.  The  latter  was  based  on  the  method  of  the  class¬ 
ical  fracture  mechanics.  Material,  geometrical  and  loading  parameters  are 
identified  explicitly  in  a  suggested  delamination  growth  model.  It  is 
therefore  tempting  to  see  whether  or  not  this  approach  remains  applicable 
to  other  types  of  delamination  processes  under  different  sets  of  material, 
geometrical  and  loading  conditions.  For  instance,  it  is  desirable  to  extend 
the  approach  to  problems  involving  compressive  and/or  fatigue  loading.  Also, 
the  growth  behavior  of  an  internal  delamination  rather  than  free  edge  delam¬ 
ination  is  another  practical  laminate  failure  problem  frequently  encountered 
in  design. 
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Thus,  it  is  the  purpose  of  the  present  study  to  expand  the  fracture 
mechanics  method  to  describe  the  mechanisms  of  delamination  in  graphite-epoxy 
laminates,  especially  under  compressive  static  and  fatigue  loads.  In  this 
effort,  four  categories  of  delamination  are  studied.  The  first  problem  is 
the  simple  case  of  free  edge  delamination  under  static  tension,  case  a.  Fig.  1. 
The  selection  of  this  problem  provides  a  ready  check  of  the  Wang/Crossman  model, 
and  helps  to  define  the  necessary  material  constants,  such  as  Gc  and  a The 
second  problem  is  the  free  edge  delamination  under  static  compressive  loading. 

In  this  case,  the  material  and  geometrical  conditions  remain  the  same  as  in 
the  first  problem;  only  the  loading  direction  is  reversed,  case  b.  The  third 
problem  is  the  delamination  emanating  from  an  implanted  interlaminar  defect.  The 
implanted  defect  is  a  thin  teflon  film  strip  (>^50  micron  thickness)  imbedded 
at  the  intended  interface  of  the  laminate,  see  case  c.  Fig.  1.  The  fourth 
problem  is  similar  except  that  the  imbedded  interlaminar  deflect  is  in  a 
circular  shape  (**  6  mm  diameter),  see  case  d.  Fig.  7. 

In  addition  to  studying  the  delamination  mechanisms  under  static  loading, 
all  four  problems  are  also  studied  under  constant  amplitude  fatigue  loadings. 

An  understanding  of  the  fatigue  failure  modes  and  the  delamination  growth 
mechanisms  is  essential  to  the  development  of  a  deterministic  cumulative  damage 
model  which  is  necessary  for  the  reliability  analysis  of  laminated  composites. 

As  far  as  delamination  mechanisms  are  concerned,  one  of  the  central  questions 
is  whether  or  not  the  mechanisms  of  damage  are  the  same  under  static  and 
fatigue  loadings.  Recent  experiments  by  Reif snider,  et  al .  [28]  on  damages 
in  laminates  indicate  that,  given  a  type  of  laminate,  there  seems  to  exist 
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a  damage  state  that  is  characteristic  to  the  particular  laminate,  whether 
it  be  under  static  or  fatigue  loading.  In  addition,  a  certain  load  threshold, 
or  a  fatigue  life  threshold  can  be  determined,  for  which  the  characteristic 
damage  state  occurs.  These  observations  suggest  that  a  law  of  deterministic 
mechanics  may  be  found,  which  relates  the  applied  load  to  the  fatigue  life  for 
a  given  state  of  sub-larainate  damage. 

In  view  of  this  background,  the  present  study  also  aims  to  use  the 
physical  parameters  identified  in  the  static  delamination  problems  in  an 
attempt  to  formulate  a  predictive  deterministic  model  for  delamination 
growth  under  fatigue  loads. 

Needless  to  say,  such  an  attempt,  though  ambitious,  is  only  exploratory 
at  this  stage.  Gains  as  well  as  open  questions  encountered  in  this  effort 
will  be  discussed  in  the  ensuing  sections. 

Section  II  discusses  the  experiment  and  the  experimental  results.  When¬ 
ever  possible,  factors  influencing  the  mechanisms  of  delamination  are  identi¬ 
fied  and  investigated  in  light  of  the  experimental  data. 

Section  III  presents  the  analytical  model  and  the  associated  numerical 
technique.  Problems  studied  in  the  experiment  are  then  simulated  numerically. 
Here,  the  limitations  of  the  present  numerical  method  for  treating  the  general 
types  of  delamination  crack  are  discussed.  Future  development  of  more 
suitable  methods  have  also  been  suggested  in  this  section. 

Finally,  a  deterministic  cumulative  fatigue  dama'ge  model  is  formulated 
in  Section  IV.  The  model  is  based  on  the  results  of  both  analysis  and 
experiment.  A  concluding  discussion  is  given  in  Section  V. 
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II.  EXPERIMENT  AND  DISCUSSIONS 
1.  Description  of  Experiments 

In  this  study  program,  the  experimental  plan  was  designed  to  include 
three  major  types  of  tests:  (a)  base-line  laminate  property  tests,  (b)  static 
load-delamination  growth  tests;  and  (c)  fatigue  load-delamination  growth 
tests. 

The  basic  material  system  chosen  for  the  experimental  plan  was  the 
Hercules  Magnamite  AS-3501-06  graphite-epoxy  system.  The  material  came 
in  prepreg  form,  and  laminates  were  fabricated  in  accordance  with  NADC  re¬ 
quirements.*  All  test  specimens  were  cut  using  a  diamond  saw,  and  end- tabbed 
with  fiber-glass  laminates.  Fig.  2  illustrates  the  general  dimensions  of  the 
tensile  specimens,  while  Fig.  3  shows  the  general  dimensions  of  the  compressive 
specimens.  In  each  case,  there  are  three  types:  one  plain  (without  Implanted 
defects)  and  two  with  implanted  defects,  as  indicated  in  Fig.  2  and  Fig.  3. 

All  tests  were  performed  on  a  closed-loop  Instron  hydraulic  testing 
machine.  A  specially  designed  test  fixture  was  used  for  all  the  compression 
tests.  Fig.  4  displays  the  schematics  of  the  compression  test  fixture.  This 
fixture  is  similar  to  that  used  in  the  Celanese  compression  test.  No  special 
fixture  was  required  for  the  tension  tests. 

The  test  environment  for  all  tests  was  maintained  at  a  room  temperature 
(  25°C)  and  ambient  humidity  ( 60%  R.H.)  condition.  Other  details  of  the 

individual  tests  are  given  in  the  following: 


*  All  laminates  were  layed  up  at  Drexel  University  and  cured  at  NADC 
facilities. 
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(a)  Base-line  Property  Tests  -  The  purposes  of  the  base-line  tests  were  to 
qualify  the  material  system  and  the  laminate  fabrication  processes,  and  to 
determine  the  material  constants  that  are  necessary  in  the  analysis  models. 

For  purpose  of  determining  the  basic  lamina  properties,  uniaxial  stress- 
strain  diagrams*  for  the  [0^2],  [90^],  [+4 ^ ] 2 s  laminates  were  obtained  under 
tension  loading,  and  for  the  [0  2],  [+45] 2S  unc*er  compressive  loading. 

Included  in  the  basic  lamina  properties  were  also  test  findings  of  [0^ 1 
fiber  contents  and  specific  gravity.  Specific  gravity  was  determined  by 
weighing  and  measuring  a  volume  of  the  material.  Fiber  content  was  deter¬ 
mined  from  photomicrographs  of  the  untested  unidirectional  material. 

Static  stress-strain  diagrams  were  also  obtained  for  the  following 
laminates;  for  static  tension. 


[02/902/452/-452]g 
[902/02/-452/452]s 
[ 0/90/0/ 90/45 /-4S/45/-45] 
[90/0/90/0/-45/45/-45/45] 
[45/— 45/0/90  ] 

L  S 

[-45/ 45/90/0  ] 

L  S 

and  for  static  compression. 


s 

s 


[02/902/452/-452]s 

[902/02/-452/452]s 

[0/90/0/90/45/-45/45/-45] 

s 

[45/-45/0/90.] 

/  s 


*  A  90°-strain  gage  rosette  was  used  to  measure  the  longitudinal  and 
transverse  strains  for  both  tension  and  compression  tests. 

8 


NADC- 79056-60 


In  all  these  static  tests,  tension  or  compression,  the  loading  was  main¬ 
tained  at  a  constant  cross-head  displacement  rate  of  0.012  cm/mln. 

Table  1  summarizes  the  individual  tests  and  the  number  of  specimens 
used  In  eacfc  test. 

(b)  Static  Load-Delamination  Growth  Tests  -  The  purpose  of  this  series  of 
tests  is  to  develop  an  understanding  of  the  load-delamination  growth  relation¬ 
ship  for  the  various  laminates  under  both  static  tension  and  static  compression 
loadings.  In  order  to  reliably  determine  the  load  at  onset  of  delamination, 
and  to  document  the  load-growth  relation,  a  step-loading  scheme  was  employed. 
Briefly,  each  specimen  is  subjected  to  step-loading  with  small  increments. 

After  each  step-load,  the  specimen  is  removed  from  the  tester  and  inspected 
by  x-radiographing  for  any  damages  that  may  occur  during  the  step-load.  The 
specimen  is  then  reloaded  for  the  next  step-loading,  x-ray  inspection,  and 
so  on. 

The  following  laminates  were  step-loaded  in  tension: 

I02/902/452/-452]a 

(902/02/-452/452]s 

(0/90/0/90/45/-45/45/-45] 

s 

[ 90/0/90/ 0/— 45/ 45/— 45 / 45] 

s 

[45/-45/0/90,] 
i  s 

[-45/45/90/0,] 

/  s 

[45/— 45/ 0/90^ 1 s  -  type  A 
[ 45/— 45/ 0/902 ] g  -  type  B 

and  the  following  laminates  were  step-loaded  in  compression: 


9 


NADC- 79056-60 


[°12] 

[02/902/452/-452]s 

[902/02/-452/452]g 

[02/902/452/-452]g  -  type  A 

[02/902/452/-4523s  -  type  C 

[0/90/0/ 90/ 45/— 45/ 45/— 45 1 s 

[ 0/ 90/0/ 90/ 45/ —45/ 45/ —4 5 ] s  -  type  A 

[0/90/0/90/ 45/-45/45/-45]g  -  type  C 

[45/-45/0/90-] 

2  s 

In  the  above,  those  laminates  designated  as  type  A  have  an  implanted 
interface  defect  consisting  of  a  0.25  in.  (6.35  mm)  Teflon  strip;  type  B 
laminates  have  a  0.5  in.  (12.7  mm)  Teflon  strip,  and  type  C  laminates  have 
a  circular  Teflon  disk  of  0.25  in.  (6.35  m»)  in  diameter.  All  defects  are 
located  at  the  mid-plane  of  the  laminates;  see  Fig.  2  and  Fig.  3. 

Table  2  summarizes  these  tests  and  the  number  of  specimens  used  in  each 

test. 

(c)  Fatigue  Life  -  De lamination  Growth  Tests  -  The  purpose  of  this  series 
of  tests  is  to  establish  the  fatigue  load-life-damage  interrelationships  of 
the  laminates  under  either  tension-tension  or  compression-compression  fatigue 
loadings.  All  tests  were  conducted  under  constant  amplitude  load  cycling  with 
a  minimum  to  maximum  load  ratio,  R  -  0.1,  at  a  frequency  of  5  Hertz. 

Again,  in  order  to  determine  the  life  (number  of  cycles)  at  onset  of  delami¬ 
nation,  and  to  document  the  life-damage-growth  history,  specimens  are 
periodically  removed  from  the  tester  and  inspected  by  x-radlography.  In  each 
type  of  laminate,  at  least  three  fatigue  load  levels  were  selected  for  the 
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fatigue  tests.  In  general,  all  the  selected  fatigue  load  levels  are  below 
the  static  load  for  which  onset  of  delamination  occurs. 

The  following  laminates  were  studied  in  tension-tension  fatigue: 

[45/— 45/ 0/ 902 ] s 

[45/— 45/ 0/ 90o ]  -  type  A 

[45/— 45/0/ 902 ] s  -  type  B 

and  the  following  laminates  were  studied  in  compression-compression  fatigue: 

[45/-45/0/90„ ] 
l  s 

t02/902M52/-452]s 

^°2^9°2^52>,-452^s  ”  type  A 

[02/902/452/-452]g  -  type  C 

[902/02/-452/452]s 

[0/90/0/90/45/-45/45/-45] 

s 

[ 0/90/ 0/90/45 /-45/45/ -45 ]g  -  type  A 
[0/90/0/90/45/-45/ 45/ —45 ] g  -  type  C 

Table  3  summarizes  these  tests,  the  number  of  fatigue  load  levels  and 
the  number  of  specimens  used  in  each  test. 

2 .  Discussion  of  Experimental  Results 

(a)  Base-line  Property  Tests  -  Results  from  all  the  static  property  char¬ 
acterization  tests  were  within  the  normal  range  of  values  for  this  material. 
Average  properties  for  each  of  the  laminates  are  summarized  in  Table  4.  In 
particular,  the  basic  properties  obtained  from  the  unidirectional  laminates 
are: 
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Tension 

Compression 

el 

140  GPa 

137  GPa 

et 

11.1  GPa 

- 

VLT 

0.29 

0.28 

VTL 

0.027 

- 

glt 

Fiber  volume,  % 

Sp.  gravity 

5.6  GPa 

66.5 

1.55 

5.4  GPa 

It  is  noted  that  the  longitudinal  Young  modulus  E  and  Poisson  ratio  v 

Li  Li 

from  both  the  tension  test  and  the  compression  test  are  essentially  the  same. 

Also,  in  obtaining  the  in-plane  shear  modulus  for  the  unidirectional  laminate, 

data  from  [+45 ]„  were  used.  Specifically,  the  average  values  for  E  and  v 

z  s  x  xy 

found  from  [+45 ] 2 s  are  use(l  in  computing  approximately  by  the  formula: 


v  ) , 

xy 


The  values  of  G  calculated  using  tension  test  results  and  using  compression 
Li 

test  results  are  also  very  close;  see  Table  4. 

However,  the  values  for  the  compressive  strength  from  Table  4  should  be 
viewed  with  caution.  They  are  uniformly  lower  than  the  tensile  strengths, 
especially  for  0°-ply  dominant  laminates.  It  is  suspected  that  the  ultimate 
compressive  laminate  failures  may  have  involved  some  degree  of  structural 
buckling. 


(b)  Static  Load-Delamination  Growth  Tests  -  In  this  series  of  tests,  eight 
different  laminates  were  tested  under  step-tension  loading,  and  nine  different 
laminates  were  tested  under  step-compression  loading.  In  each  case,  the 
various  parameters  effecting  delamination  have  been  investigated. 
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De lamination  Growth  by  Tension 

The  first  group  of  laminates  examined  were  the  quasi-isotropic  laminates, 
CO2/ 90^/ 45^/— 452 ] s  and  [90^/ 0^/— ^5^/ — 45 ] g -  In  each  case,  three  specimens 
were  tested.  The  difference  between  the  two  types  of  laminates  is  basically 
the  stacking  sequence  of  the  0°-ply  and  the  90°-ply.  Under  uniaxial  tension, 
both  laminates  suffered  transverse  cracking  in  the  90°-ply  during  the  early 
stage  of  loading.  As  the  applied  load  increased,  cracks  were  observed 
in  the  +45°  plies  at  about  80%  ultimate  load.  Edge  delamination  in  the  +45° 
interface  was  first  observable  at  about  90%  ultimate  load.  Both  laminates 
failed  at  essentially  the  same  ultimate  load  of  540  MPa  (68.4  ksi) . 

Fig.  5  displays  the  load-sequence  x-radiographs  of  a  [02/902/452/-452 ] 
laminate  under  tensile  loading.  It  is  seen  that  transverse  cracks  in  the 
90°-direction  are  visible  from  about  255  MPa  load.  The  crack  density  increases 
with  the  applied  load  until  +45°  cracks  appear  at  about  426  MPa.  These  +45° 
cracks  increased  rapidly  with  load,  causing  edge  delamination  which  initiated 
at  about  520  MPa.  The  delamination  propagated  unstably,  resulting  in  laminate 
failure  at  540  MPa. 

A  set  of  similar  load-sequence  x-radiographs  for  the  [902/02/-452/452]g 
laminate  is  shown  in  Fig.  6.  Here,  it  is  observed  that  onset  of  90°  trans¬ 
verse  cracks  occurred  at  about  234  MPa,  onset  of  +45°  cracks  occurred  at 
about  380  MPa,  and  onset  of  edge  delamination  at  442  MPa.  The  final  laminate 
failure  occurred  at  540  MPa;  this  is  the  same  load  as  in  the  [02/902/452/-452 ] 
laminate  shown  in  Fig.  5.  A  comparison  between  these  two  sets  of  results  is 
given  in  the  following: 


Onset  stress,  MPa 


laminate 

90°-crack 

+45 “-crack 

edge  delam. 

final  failure 

[02/902/452/-452]8 

255 

426 

520 

540 

[902/02/-452/452], 

234 

380 

442 

540 
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Clearly,  the  onset  stresses  for  the  different  cracking  events  are 
uniformly  lower  in  the  second  laminate  than  in  the  first.  Also,  it  seems 
that  growth  of  edge  delamination  is  the  cause  of  final  failure  for  both 
laminates.  The  difference  exhibited  by  these  two  sets  of  results  are  ex¬ 
plained  by  the  energy  analysis,  which  is  to  be  discussed  in  Section  III. 

The  second  group  of  laminates  examined  were  [ 0/ 90/ 0/ 90/45/— 45 / 45/ —45 ) s 
and  [90/0/90/0/-45/45/-45/45]g.  These  are  simply  the  same  two  quasi-isotropic 
laminates  examined  previously  with  their  stacking  plies  dispersed,  so  that  no 
two  plies  of  the  same  fiber  orientation  are  stacked  together.  By  dispersing 
the  stacking  plies,  the  laminate  stiffness  properties  do  not  change,  see 
Table  4.  However,  the  behavior  of  the  sub-laminate  cracking  events  is  changed 
by  ply-dispersion.  In  the  step-tension  tests,  both  laminates  showed  90°- 
cracking  at  the  applied  load  of  434  MPa;  neither  +45°-cracking  nor  edge 
delamination  was  observed  before  the  laminate  final  failure.  Furthermore,  the 
final  laminate  failure  load  for  both  laminates  was  over  620  MPa,  or  15%  higher 
than  the  undispersed  laminates. 

These  results  show  clearly  the  influence  of  ply  thickness  on  the  sub¬ 
laminate  cracking  behavior  as  well  as  the  final  strength  of  the  laminate. 
Energy  analysis  on  this  subject  is  also  included  in  Section  III. 

The  next  group  of  laminates  tested  were  the  [+45/0/ 90^ ] g  series.  Three 
types  were  examined:  the  plain  (without  implant),  the  type  A  (with  1/4"  Teflon 
strip  implanted)  and  the  type  B  (with  1/2"  Teflon  strip  implanted);  see 
dimensions  given  in  Fig.  2. 

When  loaded  in  tension,  tensile  edge  stress  is  induced  throughout  the 
90°-plies.  This  a stress  is  thus  responsible  for  edge  delamination  in  the 
90°-plies.  Test  results  show  the  following  cracking  activities: 
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90°-crackings  occurred  at  about  131  MPa  in  the  plain  laminates, 
and  occurred  at  about  103  MPa  in  both  type  A  and  type  B  laminates.  Edge 
delamination  became  observable  at  about  223  MPa  in  the  plain  laminates,  and 
about  180  MPa  in  type  A  and  type  B  laminates.  All  three  types  had  the  same 
average  failure  stress  of  about  438  MPa.  These  results  are  summarized  in  the 
following 

_ onset  stress,  MPa _ 

type  o  animate  90°-cracking  edge  delam.  final  failure 


[+45/0/90] 

—  s 

131 

223 

443 

type  A 

103 

179 

436 

type  B 

103 

172 

435 

Clearly,  the  presence  of  interlaminar  defects  in  the  type  A  and  type  B 
laminates  has  some  measurable  adverse  effect  on  the  initiation  of  transverse 
cracking  and  edge  delamination.  But  the  final  failure  load  of  the  laminates 
is  not  affected.  This  suggests  that  the  growth  behavior  of  the  sub-laminate 
cracking  events  is  probably  quite  stable,  and  that  the  extent  of  sub-laminate 
damage  prior  to  final  failure  is  also  the  same  in  all  three  types  of  laminates. 
Indeed,  from  the  load-sequence  x-radiographs  shown  in  Figs.  7,  8  and  9,  it  is 
seen  that  in  all  cases,  free  edge  delamination  growth  is  extremely  stable. 

The  influence  of  the  implanted  defect  becomes  insignificant  as  the  applied 
load  approaches  its  ultimate. 

Fig.  10  illustrates  graphically  the  90°-cracking  and  edge  delamination 
growth  as  a  function  of  the  applied  load.  It  is  seen  that  the  planted  inter¬ 
laminar  flaws  had  little  effect  on  the  growth  behavior.  It  is  essentially 
negligible  as  the  applied  load  reaches  the  ultimate. 

Again,  energy  analysis  of  the  cracking  events  for  this  series  of 
laminates  is  included  in  Section  III. 
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Finally,  several  specimens  of  [+45/90/02]s  configuration  were  also 
tested.  An  interlaminar  edge  stress  analysis  would  show  that  the  edge  normal 
stress  a is  compressive.  Thus,  no  delamination  is  expected  in  such  a 
laminate  configuration.  In  the  step-loading  experiment,  some  90°-cracking  was 
observed  from  about  820  MPa  which  is  88%  the  ultimate  strength  of  the  laminate 
(o  =  935  MPa).  Even  prior  to  final  failure,  there  was  no  noticeable  +45°- 
cracking  or  free  edge  delamination  observed. 

Table  5  summarizes  all  the  test  results  obtained  for  the  eight  types  of 
laminates.  Major  findings  from  this  series  of  tests  may  be  briefly  stated 
as  follows: 

(1)  Between  the  laminates  ^s  an<*  f^02/0 2/~45  ^ ^2^  s’  very 

little  difference  exists  in  terms  of  laminate  behavior,  including  their 
stiffnesses  and  the  stress  states.  Initiation  of  90°-cracking 
occurred  early  in  the  loading  for  the  second  laminate.  This  sub¬ 
laminate  damage  event  influenced,  in  turn,  an  early  onset  of  +45°- 
cracking  and  edge  delamination  in  the  +45°  interface.  Although  the 
ultimate  load  for  both  laminates  was  the  same,  the  process  of  damage 
development  before  final  failure  is  not. 

(2)  Between  the  laminates  [0o/90„/45./-45_ ]  and  [0/90/0/90/45/-45/45/-45]  , 

2  2  2  2  S  S 

the  difference  exhibited  in  sub-laminate  damage  development  is  significant 

The  following  comparison  shows  their  strikingly  different  behavior 


laminate 


_ Onset  stress,  MPa _ 

90°-cracking  +45°  cracking  Edge  Delara. 


Final  failure 


[02/902/452/-452]s 


[0/90/0/90 
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The  effect  of  ply  thickness  on  both  crack  initiation  and  growth  behavior 
is  clearly  demonstrated  here.  It  is  an  important  parameter  in  the 
fracture  mechanisms  associated  with  these  cracking  events. 

(3)  For  the  [+45/0/902^  series,  the  effect  of  implanted  interlaminar  defect 
is  relatively  unimportant  if  the  laminates  are  under  tension  (static 
loading).  The  defects  do  cause  early  initiation  of  damage  growth, 
emanating  from  the  defect.  This  could  be  a  factor  in  situations  where 
unstable  damage  growth  is  dominant,  however. 

Delamination  Growth  by  Compression 

Two  major  groups  of  laminates  were  tested  by  static  step-compression. 

The  first  group  consisted  of  the  [02/,902^52^-452^s  lanlinates»  and  included 

those  with  a  Teflon  strip  implant  (type  A),  and  those  with  a  Teflon  circular 

disk  (type  C)  implanted  in  the  raid-plane.  The  purpose  of  testing  this  grouped- 

♦ 

ply  laminates  was  simply  to  see  the  effect  of  the  implanted  defects  on  ■* 

delamination  under  compressive  loading. 

The  second  group  of  laminates  were  the  "dispersed"  laminates;  namely,  the 

[0/90/0/90/45/-45/45/-45]g  and  included  both  type  A  and  type  C  implants.  Thus, 

between  the  two  groups,  the  ply-grouping  effect  on  delamination  can  be 

evaluated  for  the  compression  loading  cases. 

In  addition,  some  [0lo],  [90o/0  /-450/450]  and  [+45/0/90-]  laminates 

14  Li  Li  s  —  is 

were  also  tested  by  step-compression.  For  [O^l*  was  hoped  to  detect  any 
fiber  splitting  during  the  step-loading,  or  any  other  matrix-dominated  crack¬ 
ing  prior  to  final  failure.  As  for  the  [902/02/-452/452]g,  its  behavior  should 
be  identical  with  the  [ 0^/902 ^^52/~ ^5^^  under  compression.  And,  of  course, 
the  [+45/0/902^  should  exhibit  no  delamination  under  such  loading,  as  well. 
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Figs.  11,  12  and  13  show,  respectively,  the  load-sequence  x-radiographs 
of  the  [Q ^1 $0 ^ ^ 2^ s*  type  A  and  type  C  laminates.  Free  edge  delamination 
was  observed  at  various  loading  levels,  depending  on  the  type  of  laminates. 

The  following  shows  the  onset  loads  at  delamination  and  failure: 


laminate 

[02/902/452/-452]s 
[02/902/452/-452Js  -  A 
[02/902/452/-452]g  -  C 

[902/02/452/-452]s 


Onset  stress,  MPa 
Edge  delam.  final  failure 

339.2  369.5 


131.0  286.0 


252.0  331.0 


316.0  320.0 


From  the  above,  it  is  seen  that  in  the  plain  [0o/90o/45_/-45„ ]  and 

lit  is 

[90„/0./45o/-45_]  laminates,  onset  of  edge  delamination  led  immediately  to 
1  l  l  l  s 

final  laminate  failure;  the  growth  of  delamination  was  unstable. 

In  the  type  A  laminates,  the  implanted  teflon  strip  is  exposed  at  the 
free  edges*  This  defect  caused  an  early  but  stable  edge  delamination  growth, 
as  is  seen  in  the  x-ray  pictures  depicted  in  Fig.  12.  This  growth  also  pre¬ 
cipitated  a  lowered  final  failure  stress  for  the  laminates  (about  25%  lower). 

The  implant  in  the  type  C  laminate  is  in  the  interior;  it  has  no  effect 
on  the  free  edges  of  the  laminates.  Initially,  there  are  two  independent  delam¬ 
ination  initiators,  the  implant  and  the  free  edges.  Apparently,  the  implanted 
circular  disk  had  less  effect  on  delamination  than  the  free  edge.  As  is 
seen  in  the  load-sequence  x-ray  pictures  in  Fig.  13,  edge  delamination  was 
the  first  event  to  occur.  Once  edge  delamination  started,  however,  the 
internal  circular  patch  affected  the  growth  behavior,  resulting  in  a  lowered 
final  failure  load  for  the  laminates. 
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The  delamination  growth  behavior  of  these  three  types  of  laminates  is 
illustrated  in  Fig.  14.  where  the  amount  of  delamination  growth  (in  %  of 
delaminated  area)  is  plotted  against  the  applied  load.  It  is  seen  that 
delamination  growth  in  the  ^-aminates  almost  catastropic 

following  the  initiation.  Type  C  laminates  had  a  relatively  stable  growth, 
while  type  A  laminates  had  the  most  stable  growth.  The  effect  of  the  inter¬ 
laminar  defect  on  delamination  was  significant.  This  is  in  contrast  with 
the  results  from  tensile  step-loading,  where  the  effect  of  interlaminar 
defects  on  delamination  was  relatively  insignificant;  see  Fig.  10. 

Figs.  15,  16  and  17  display  the  load-sequence  x-radiographs  for  the  three 
types  of  laminates  in  the  second  group.  Recall  that  in  this  group,  the  ply 
stacking  sequence  of  the  previous  laminate  is  dispersed.  The  following  table 
shows  the  comparisons  of  the  onset  loads  of  delamination  and  final  failure  for 
these  laminates: 


laminates 

Onset 

stress,  MPa 

edge  delam. 

final  failure 

[0/90/0/90/ 45/-45/45/-45] 

s 

379.2 

425.4 

[0/90/0/90/45/-45/45/-45]s  -  A 

138.0 

382.6 

[ 0/90/0 /90/45/-45/ 45/-45 ]g  -  C 

265.4 

365.4 

Again,  it  is  seen  that  edge  delamination  was  the  initial  damage  event  pre¬ 
cipitating  failure  in  all  the  laminates.  The  effect  of  an  interlaminar  defect  on 
the  initiation  of  delamination  remains  as  significant  as  in  the  previous  un¬ 
dispersed  laminates.  Fig.  18  illustrates  the  various  delamination  growth 

developments  as  functions  of  the  applied  load. 

Between  the  two  groups  of  laminates,  the  undispersed  and  the  dispersed, 

the  important  difference  is  again  caused  by  the  parameter  of  ply  thickness. 

Fig.  19  shows  the  delamination  growth  behavior  between  the  undispersed  and 

the  dispersed  laminates. 


19 


NADC- 79056-60 


Table  6  summarizes  all  the  results  obtained  In  this  series  of 
tests.  It  is  clear  that  interlaminar  defects  have  a  significant  adverse 
effect  on  delamination  when  the  laminates  are  subjected  to  compression.  In 
addition,  the  ply  thickness  parameter  continues  to  play  an  important  role 
in  the  delamination  growth  processes . 

(c)  Fatigue  Life  -  Delamination  Growth  Tests  -  In  this  series  of  tests,  both 
tens ion- tension  fatigue  and  compression-compression  fatigue  tests  were 
conducted.  Details  of  test  results  are  discussed  in  the  following. 


Tension-Tension  Fatigue  Tests 


In  the  tens ion- tens ion  fatigue  experiment,  only  the  [+45/0/902^  and 

the  associated  type  A  and  type  B  laminates  were  tested.  For  the  [+4 5/0/90^^ 

laminates,  five  fatigue  load  levels  were  selected.  These  were  one  load  level 

below,  three  load  levels  above,  and  one  load  level  at  about  the  static 

onset  load  of  transverse  cracking  in  the  90°-plies.  The  purpose  of  these 

tests  was  to  establish  some  quantitative  damage  growth  rates  in  terms  of  the 

fatigue  load  and  the  fatigue  life.  A  typical  life-sequence  of  x-radiographs 

for  the  [+45/0/90- ]  laminate  tested  under  ae  =  189  MPa  is  shown  in  Fig.  20. 

—  i  s  i 

Since  the  maximum  fatigue  load  is  higher  than  the  static  onset  load  for 
transverse  cracking  (at  131  MPa),  it  is  seen  from  the  x-radiographs  that 
multiple  90°- cracks  occur  very  early  during  the  fatigue  loading.  As  the 
fatigue  cycling  increases,  the  crack  density  also  increases.  At  about 
1000  cycles,  free-edge  delamination  becomes  identifiable.  The  delamination 

then  grows  with  increasing  fatigue  cycles. 

Fig.  21  shows  the  growth  of  the  90°-cracks  as  a  function  of  fatigue 


cycles,  while  Fig.  22  displays  the  delamination  growth.  It  is  seen  from  these 
figures  that  free  edge  delamination  at  all  fatigue  load  levels  occurs  whenever 
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the  transverse  crack  density  reaches  8  '“lO  cracks  per  cm.  From  there  on,  edge 
delamination  becomes  the  dominant  damage  mode,  which  leads  steadily  to  final 
laminate  failure.  Both  Fig.  21  and  Fig.  22  suggest  a  certain  damage 

growth  rate  dependent  upon  the  maximum  fatigue  load,  a^,  and  the  cyclic  loading 
history,  N.  The  exact  form  for  such  a  growth  model  will  be  discussed  in 
Section  IV. 

The  effect  of  interlaminar  defects  on  this  laminate  series  under  tension 
fatigue  may  be  determined  from  the x-radiographs  taken  for  a  [+1*5/0/90^]  -  A 
laminate  and  a  (+45/O/9O2]  -  B  laminate,  shown  respectively  by  Fig.  23  and 
Fig.  24.  Both  laminates  were  cycled  under  a  maximum  fatigue  load  of  189  MPa. 

It  is  observed  from  both  cases  that  90°-cracking  occurs  early  in  the  fatigue 
life,  followed  promptly  by  edge  delamination  eminating  from  the  inter¬ 
laminar  flaw.  Thus,  the  presence  of  different  interlaminar  flaws  had  only 
a  small  effect  on  the  growth  of  90°-cracking,  but  precipitated  early  delami¬ 
nation  initiation.  Fig.  25  shows  the  delamination  growth  under  the  same 
fatigue  load  (189  MPa)  for  the  three  types  of  laminates.  Clearly,  the  effect 
of  implanted  flaws  is  significant  enough  that  it  should  not  be  disregarded, 
although  its  effect  on  static  tension  strength  was  found  to  be  minimal. 

Results  obtained  from  tests  at  other  fatigue  load  levels  are  summarized 
in  Table  7. 

Compression-Compression  Fatigue  Tests 

As  in  the  case  of  the  static  compression  tests,  two  major  laminate 

series  were  also  tested  under  compression-compression  fatigue  loads.  These 

two  series  were  (a)  the  f 0^/902/ 45^/— 452 ] s  and  the  associated  type  A  and 

type  C  laminates,  and  (b)  the  dispersed  [0/90/0/90/ 45/-45/45/-45J  and  the 

s 
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associated  type  A  and  type  C  laminates.  In  addition,  laminates  of  [+45/0/902 J 
and  [9O2/O2/-452/452 ]  were  also  tested.  In  each  case,  two  or  three  fatigue 
load  levels  were  selected;  see  Table  8. 

For  the  [0^ /902/^52 /— A5^ ] s »  type  A  and  type  C  series,  typical  x-radio- 
graphs  are  displayed,  respectively,  in  Figs.  26,  27  and  28.  These  were  all 
under  a  222  MPa  maximum  fatigue  load.  This  load  is  much  higher  than  the 
static  onset  load  for  delamination  for  type  A  laminate  (131  MPa) ,  but  is  lower 
than  that  for  the  plain  and  the  type  C  laminates  (339  MPa  and  252  MPa,  re¬ 
spectively)  .  Thus,  in  the  case  of  type  A  laminate,  edge  delamination  was 
found  during  the  first  cycle  of  loading.  Growth  of  delamination  was  rather 
stable  until  fatigue  life  exceeds  10^  cycles.  For  the  plain  and  the  type  C 
laminate,  delamination  occurred  later  in  the  fatigue  life  but  with  a  faster 
growth  rate. 

Fig.  29  shows  the  delamination  growth  curves  obtained  from  three  fatigue 
load  levels  for  the  plain  laminates.  The  growth  rates  in  these  curves  are 
higher  compared  with  those  found  under  tension  fatigue.  The  same  growth 
behavior  is  observed  also  for  the  type  A  and  type  C  laminates,  see  Fig.  30 
and  Fig.  31,  respectively.  Both  show  an  accelerated  growth  rate  under  com¬ 
pression  fatigue.  This  result  is  thought  to  be  due  to  the  buckling  action  of 
the  test  specimen  under  compression.  This  behavior  must  be  considered  in  any 
growth  model  formulation. 

The  effect  of  implanted  interlaminar  flaws  is  illustrated  in  Fig.  32. 
Under  a  221  MPa  fatigue  load,  it  is  seen  that  edge  delamlnation  initiates 
at  different  lives  depending  on  the  type  of  flaw  implanted.  For  instance, 
type  A  laminate  suffers  delamination  at  the  first  cycle,  type  C  begins  to  show 
delamlnation  at  about  200  cycles,  while  the  plain  laminate  shows  delamination 
at  a  life  exceeding  2000  cycles.  The  rate  of  growth  in  all  three  cases  takes 


22 


NADC- 79056-60 


a  sharp  increase  at  about  25%  delamination,  and  the  fatigue  lives  of  all 
three  cases  remain  within  one  decade  of  variation  of  each  other. 

The  purpose  of  testing  the  [0/90/0/90/45/— 45/ 45/— 45 ]  and  the  associated 

s 

type  A  and  type  C  laminates  was  again  to  assess  the  effect  of  ply-grouping 
when  the  loading  condition  is  compression-compression  fatigue.  In  the  exper¬ 
iments  on  this  series  of  laminates,  the  same  general  behavior  was  found  as  in 
the  previous  series.  In  all  cases,  edge  delamination  was  induced  under  the 
fatigue  loads,  and  while  the  effect  of  the  implanted  flaw  was  noticeable,  it 
was  not  significant. 

Figs.  33,  34,  and  35  illustrate  respectively,  the  growth  of  delamination 
by  a  life-sequence  of  x- radiographs  for  the  plain  (under  276  MPa), 
type  A  and  type  C  (both  under  219  MPa)  laminates.  From  these  x- radiographs, 
delamination  growth  curves  can  be  obtained.  Fig.  36  shows  a  typical  set  of 
delamination  growth  curves  for  the  type  A  laminate  under  three  different 
fatigue  load  levels.  It  is  noted  that  these  curves  are  similar  with  those 
found  earlier  for  the  [0o/90„/45o/-45o ]  series,  except  that  delamination 

4  L  L  i  s 

growth  is  suppressed  by  the  dispersion  of  the  plies. 

In  order  to  see  the  ply-grouping  effect  more  quantitatively,  examine 
the  comparative  growth  curves  shown  in  Figs.  37,  38,  and  39.  In  all 
cases,  the  dispersed  laminates,  whether  they  contain  implanted  flaws  or  not, 
display  a  much  "stronger"  behavior  against  delamination  than  their  undispersed 
counterparts.  The  dispersed  laminates  show  more  than  2  decade  longer  life 
than  the  undispersed.  This  significant  simple  geometrical  parameter  is  seen 
to  have  a  dominant  effect  on  the  failure  mechanisms  of  the  laminates. 
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Finally,  teat  results  for  the  [+45/0/90^ ls  and  [9O2/O2/-452/452 J8 
laminates  were  as  expected.  The  [+45/0/902 3 s  failed  in  fatigue  without 
any  edge  delamination,  and  the  Shaved  similarly  to 
[02/902/*52/-»5jJ, 

All  these  results  are  summarized  in  Table  8. 
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HI.  AN  ENERGY  METHOD  FOR  SUB-LAMINATE  CRACK  GROWTH 

In  this  section,  the  approach  being  taken  in  the  development  of  a  crack 
growth  model  will  be  described.  To  this  end,  a  brief  review  is  presented 
describing  the  major  sub-laminate  crack  mode  and  the  classical  fracture  con¬ 
cept  of  strain  energy  release  rate.  Finally,  the  cracking  events  documented 
from  experiments  will  be  analyzed  using  the  energy  model  presented  in  this 
section. 

1 .  Physical  Processes  of  Transverse  Cracking  and  Edge  Delamination 
Examination  of  the  earlier  test  results  show  that  matrix  and/or 
matrix-fiber  interface  damage  almost  always  precedes  the  laminate  failure. 
Particularly,  90°-cracking  (transverse  cracks)  and  free-edge  delamination  have 
been  identified  as  the  two  major  damage  modes  which  precede  the  final  failure 
in  the  laminates  studied  in  this  report.  In  order  to  properly  model  each  of 
these  cracking  processes,  a  closer  analysis  of  the  mechanisms  of  these  events 
becomes  essential. 

(a)  Transverse  Cracking  Process  Transverse  cracks  form  in  plies  when  a 
sufficiently  large  tensile  stress  exists  transverse  to  the  fibers.  Consider, 
for  example,  a  cross-ply  (0/90)g  laminate  which  is  subjected  to  tension  in  the 
0°-direction;  at  some  critical  stress  level,  the  90°-ply  begins  to  suffer 
transverse  cracks.  Laboratory  experience  indicates  the  following  characteris¬ 
tics  of  the  transverse  cracking  process: 

Reproducible  Onset  Load  In  a  given  population  of  identical  specimens  sub¬ 
jected  to  identical  test  conditions,  the  critical  load  at  initiation  of 
transverse  cracking  is  surprisingly  reproducible,  especially  in  view  of  the 

large  scatter  experienced  in  material  properties  (such  as  the  transverse  strength 
of  a  90 laminate ) .  This  suggests  that  the  transverse  crack  initiation 

mechanisms  are  predominantly  influenced  by  the  geometrical  conditions  surround¬ 
ing  the  crack,  rather  than  by  the  material  property  variations. 


NADC- 79056-60 


Thickness  Effect  and  Shear  Lag  X-radiographs  in  earlier  experiments  showed 
that  transverse  cracks  in  the  90°-ply  generally  follow  the  matrix-fiber  inter¬ 
face;  in  a  trans-thickness  view,  such  as  sketched  in  Fig.  40,  the  crack  is 
believed  to  form  inside  the  90°-ply,  and  extend  vertically  across  the  90°-ply, 
eventually  terminating  at  the  0°-90°  interface.  Clearly,  the  outside  0°-plies 
must  have  a  restraining  and  arresting  effect  as  the  transverse  crack  approaches 
this  interface. 

If  strain  energy  is  released  as  the  transverse  crack  propagates,  it  is 
easy  to  see  that  the  farther  a  transverse  crack  is  allowed  to  propagate  (i.e., 
the  thickness  of  90°-ply  is  larger)  the  greater  the  strain  energy  released 
by  the  propagating  crack.  Thus,  from  the  energy  point  of  view,  a  thin 
90°-ply  layer  can  provide  only  a  small  amount  of  energy  release.  Consequently, 
a  higher  applied  load  is  required  to  propagate  a  crack  in  the  thin  90°-ply 
layer.  This,  of  course,  is  consistent  with  the  experimentally  observed 
thickness  effect,  discussed  extensively  in  Section  II. 

After  the  transverse  crack  is  arrested  at  the  0°/90°  interface,  the  load 
normally  carried  by  the  90®-ply  before  the  crack  must  now  be  transmitted 
through  the  0°/90°  interface,  creating  a  local  region  of  interlaminar  shear 
stress  distribution,  called  the  shear  lag.  Obviously,  the  size  of  the  shear 
lag  depends  again  on  the  thickness  of  the  90°-ply  layer.  Generally,  it  is 
about  4-times  the  thickness  of  the  90°-layer. 

Theoretically,  if  the  material  properties  of  the  90°-layer  are  perfectly 
uniform,  then  any  two  adjacent  transverse  cracks  must  maintain  a  minimum  dis¬ 
tance  equal  to  the  shear  lag.  Thus,  at  the  onset  of  transverse  cracking,  one 
would  actually  see  multiple  cracks  that  are  separated  by  the  shear-lag  spacing. 
This  has  been  referred  to  as  the  "characteristic  spacing"  by  Reif snider. 
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In  actuality. however,  the  variability  of  the  material  properties  distortes 
the  so-called  "characteristic  spacing."  The  growth  of  the  transverse 
cracking  is  a  continuous  function  of  the  applied  load,  as  observed  in  the 
earlier  experiments.  The  lack  of  a  "characteristic  spacing"  complicates  the 
definition  of  cracking  initiation;  and  necessitates  a  statistical  approach 
to  address  the  random  nature  of  the  material  variability. 

Effect  of  Transverse  Cracking  on  Edge  Delamination  The  aforementioned  shear- 
lag  region  may  suffer  local  damage  in  the  form  of  delamination,  due  to  the 
concentration  of  interface  shear  stress  near  the  transverse  crack-tip.  Since 
the  size  of  the  shear-lag  is  proportional  to  the  90°-ply  layer  thickness,  the 
likelihood  of  transverse  cracking  induced  delamination  also  increases  with  the 
90°-layer  thickness.  Thus,  a  laminate  which  is  already  delamination  prone  due 
to  edge  effects  will  suffer  a  compounded  delamination  growth  if  it  under¬ 
goes  simultaneously  transverse  cracking.  Such  an  interacting  effect  is  3- 
dimensional  in  nature,  and  requires  3-dimensional  crack -mode ling. 

(b)  Edge  Delamination  Process  As  has  been  discussed  earlier  in  Section  I, 

edge  delamination  stems  from  the  boundary  edge  stresses  which  exist  along  ply 

interfaces  near  the  laminate  free  edges.  These  interlaminar  stresses 

can  be  computed  accurately  by  a  stress  analysis  based  on  ply  elasticity. 

Again,  consider  as  an  example  the  (0°/90°)g  laminate  subjected  to  uniaxial 

tension.  An  edge  stress  analysis  will  show  (see  e.g.,  Ref.  14)  that  a  tensile 

az  stress  exists  throughout  the  90°-ply  layer,  and  a  combination  of  tensile  a z 

and  shearing  t  exists  along  the  0°/90°  interface.  Under  these  stresses,  a 
xz 

mode-I  crack  along  the  mid-plane  (90°/90°  interface),  or  a  mixed-mode  crack 

along  the  0°/90°  interface  are  both  possible,  depending  on  which  of  the  two 

possible  sites  meet  the  cracking  conditions. 
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Furthermore,  the  edge  stresses  exist  only  within  a  small  boundary  edge  region. 
The  size  of  the  boundary  region  is  again  found  to  depend  on  the  thickness  of  the 
layers  along  whose  interface  the  edge  stresses  are  distributed.  Hence,  the 
thicker  the  material  layers,  the  larger  the  size  of  the  boundary  region.  If 
an  interlaminar  crack  is  to  form  in  this  region,  it  may  be  inferred  again  that 
the  available  strain  energy  released  during  the  formation  of  the  crack  is 
limited  by  the  layer  thickness.  Or,  alternatively,  the  critical  load  at  onset 
of  edge  delamination  depends  on  the  material  layer  thickness .  This  is 
also  clearly  consistent  with  the  experimental  results. 

Once  the  free  edge  crack  is  formed,  the  boundary  region  propagates  inward 
to  the  interior  of  the  laminate.  If  the  width  of  laminate  is  large  compared 
to  the  size  of  crack,  stable  growth  will  follow.  Otherwise,  the  growth  may 
become  unstable. 

The  foregoing  physical  analyses  are  based  primarily  on  experimental 
observations.  With  this  understanding,  the  development  of  a  descriptive  model 
will  be  presented  next,  on  the  basis  of  the  energy  release  rate  concept. 

2 .  The  Strain-Energy  Release  Rate  Concept 

The  concept  of  the  strain  energy  release  rate  as  a  fracture  growth  criterion 
for  brittle  cracks  in  elastic  solids  is  of  long  standing.  In  his  first 

paper  published  in  1920,  Griffith  [29]  attributed  the  ultimate  tensile  strength 
of  glass  to  the  propagation  of  crack-like  surface  flaws;  stress  concentration 
surrounding  the  worst  (largest)  flaw  cause  unstable  propagation.  He 
postulated  that  a  change  (increase)  in  the  crack  size  causes  a  decrease  in  the 
stored  strain  energy  near  the  crack-tip;  at  the  same  time,  the  solid  body 
increases  its  surface  area  so  that  the  total  surface  energy  is  also  increased. 
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Let  U  be  the  strain  energy  and  S  the  surface  energy;  then  for  an  initial  flaw 
of  size  a,  a  crack  growth  is  induced  if 


3U  JS 
3a  ^  3a 


(3-1) 


The  quantity  3U/3a  depends  on  the  stress  field  near  the  crack-tip  which 
in  turn  depends  on  the  geometry  of  the  crack  and  the  magnitude  of  the  applied 
load.  This  quantity  represents  the  driving  force,  and  is  commonly  referred  to 
as  the  available  energy  release  rate,  G(a).  The  quantity  3S/3a  depends  on 
the  microstructure  of  the  material  and  the  magnitude  of  the  molecular  bond; 
it  is  thus  considered  a  material  property,  representing  the  resistance  of 
the  material  against  crack  growth.  It  is  also  commonly  called  the  critical 
energy  release  rate,  G^.  For  most  materials,  is  a  constant,  independent 
of  the  crack  size  or  the  applied  load  (but  may  depend  on  the  crack  direction) . 

The  Griffith  crack  growth  criterion,  (3-1)  may  be  rewritten  in  the  form 


G  >  G  . 
—  c 


(3-2) 


(a)  The  Crack-Closure  Representation  From  an  elastic  energy  balance  con¬ 
sideration,  Irwin  [25]  showed  that  the  strain  energy  released  in  an  incre¬ 
mental  crack  extension  can  be  equated  to  the  amount  of  work  done  required  to 
close  it  again.  Let  Aa  represent  an  infinitesimal  crack  extension  from  the 
initial  crack  of  size  a,  as  shown  in  Fig.  4 1.  Let  Au  be  the  relative  dis¬ 
placements  between  the  mating  surfaces  along  Aa  when  the  crack  extension  is 
open,  as  in  Fig.  41(b);  let  a  be  the  surface  stresses  along  Aa  when  the  crack 
extension  is  closed,  as  in  Fig.  41(c) 

Then,  the  work  done  in  closing  the  crack  extension  Aa  can  be  expressed  as 

AW  =  —  I  o’* Ail  da  (3-3) 

1  J0 
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By  Irwin's  elastic  energy  balance,  the  strain  energy  release  rate  correspond¬ 
ing  to  the  crack  of  initial  size  a  is  given  by 


G 


lim  1  (Aa 
Aa-»0  2Aa  j  ^ 


M  % 

a*Au  da 


(3-4) 


Substituting  the  normal  and  shearing  stress  components  of  the  surface 
stress  a,  and  the  relative  displacement  Au  into  (3-4),  the  modal  components 
of  G  are  obtained  as 


lim  1 
Aa->0  2Aa 


A  a 


o  Av  da 


=,  11m  1 

’ll  =  Aa-+0  2Aa 


„  lim  1 

’ill  Aa-*)  2Aa 


•A  a 

0 

Aa 


t  Au  da 
xy 


(3-5) 


yz 


Aw  da 


where  G^,  G^  ^  and  are,  respectively,  the  energy  release  rates  correspond¬ 

ing  to  the  mode-I,  mode-II  and  mode- III  crack  opening,  such  as  illustrated  in 
Fig.  42.  The  total  energy  release  rate  of  the  generally  mixed-mode  crack  ex¬ 
tension  is  given  by 


GI  + 


GII  +  GIII 


(3-6) 


(b)  The  Finite-Element  Scheme  In  order  to  evaluate  the  integrals  in 
(3-4)  or  (3-5),  the  exact  stress  field  near  the  crack-tip  must  be  determined. 
But,  the  stress  field  surrounding  a  sharp  crack  is  singular  in  nature;  and 
in  a  generally  anisotropic  solid,  the  mode-I,  mode-II  and  mode-III  cracking 
actions  are  mutually  coupled.  Thus,  the  crack-closure  representation  of 
Irwin's  must  be  adapted  by  applying  a  finite  element  approximation  for 
problems  involving  complicated  layering  structures.  A  technique  of 
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nodal-release-and-close  has  been  implemented  recently  by  Wang  and  Crossraan 
[20]  to  approximate  the  Irwin  integrals  of  (3-5).  Fig.  43  illustrates  the 
finite  element  representation  of  the  crack-tip  region  shown  in  Fig.  41(a). 

Here,  the  crack  of  size  a  is  shown  with  the  crack-tip  at  node  C.  Under  a 
given  far  field  loading  condition,  the  displacements  of  node  C,  u^,  can  be 
determined  by  the  finite  element  solution.  Then,  a  small  incremental  crack 
extension  (self-similar,  for  instance)  Aa  is  introduced  by  opening  up  node  C, 
with  node  C  becoming  two  separate  nodes  f  and  g,  as  shown  in  Fig.  43(b).  A 
new  round  of  finite  element  solutions  is  obtained,  given  the  new  crack  geometry 
(the  far  field  loading  condition  remains  unchanged) .  This  provides  the  nodal 
displacements  of  node  f  and  node  g,  denoted  respectively  by  u^  and  u^.  The 
last  step  is  to  force-close  node  f  and  g  to  return  to  the  original  position 
of  node  C  by  prescribing  the  common  displacement  uc  for  both  node  f  and  node  g. 
The  finite  element  solution  of  the  last  step  gives  the  necessary  nodal  force 

rv 

Fc  required  for  the  closure.  The  work  done  in  the  process  of  a  finite  crack- 
closure  is  given  by 


AW  *  y  [Fx(uf"ug)  +  Fy(vf"vg)  +  FZ(wf'wg)] 


(3-7) 


Thus,  for  approximation,  the  strain  energy  release  rate  G  may  be 
expressed  as  * 


G  =  AW/Aa. 


(3-8) 


*  The  adequacy  of  this  kind  of  approximation  has  been  a  source  of  great  con¬ 
cern.  A  more  detailed  assessment  of  such  approximation  is  found  in  Ref.  [20]. 
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Consequently,  the  ntodal  components  of  G  corresponding  to  the  Initial  crack 
of  size  a  are  represented  by 

G_  ■  F  (v  -v  )/2da 

I  y  y  g 

Gn  "  Fx<uf“ug)/2Aa  (3-9) 

GIII  *  Fz(wf-wg)/2Aa 

By  successively  opening  and  closing  the  nodes  along  a  prescribed  crack 
path,  the  energy  release  rate  G  can  be  computed  as  a  function  of  the  extend¬ 
ing  crack's  size,  i.e.,  G  =  G(a)  under  a  given  far-field  loading  condition. 

(c)  Sub-laminate  Crack  Initiation  and  Growth  Criteria  In  this  study,  the 
only  sub-laminate  cracks  investigated  are  the  90°-cracking  and  the  free-edge 
delamination.  In  order  to  set  up  a  crack  initiation  and  growth  criterion  for 
each  of  the  two  types  of  cracking,  it  is  assumed  that  micro-flaws  exist 
along  the  fiber-matrix  interface,  and  also  along  the  ply-to-ply  interface. 
These  flaws  are  generally  of  a  size  in  the  order  of  the  fiber  diameter.  This 
dimension  is  in  the  realm  of  micro-mechanics,  whose  crack  propagation  mech¬ 
anisms  will  not  be  dealt  with  here.  Instead,  it  is  postulated  only  that  these 
micro-flaws  would  propagate  and  coalesce  into  a  macroscopic  proportion  under 
a  certain  critical  far  field  loading  condition.  The  resulting  crack  is 
recognizable  only  at  the  macro-scale,  where  the  basic  axioms  of  ply-elasticity 
are  valid*.  The  minimum  crack  size  first  formed  at  the  macro-scale  depends  on 
the  micro-structure  of  the  material  (perhaps  also  on  a  host  of  other  material 
processing  factors,  such  as  the  curing  process,  post-cure  handling,  etc). 


*  One  of  the  basic  axioms  in  ply-elasticity  is  that  within  each  layer  the 
material  is  represented  by  a  homogeneous,  anisotropic  medium.  Material  dis¬ 
continuity  exists  only  across  the  interface  of  two  different  material  layers. 
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This  minimum  macro-crack  size  is  also  called  the  critical  flaw  size,  denoted 
by  ac.  The  value  of  ac  is  generally  unknown  and  can  be  random  in  nature.  It 
is  then  regarded  as  a  material  property  at  the  macroscopic  level. 

The  Griffith  criterion  of  (3-2),  when  applied  to  predict  the  initiation 
of  a  sub-laminate  cracking,  is  replaced  by 

G(a  )  >  G  .  (3-10) 

c  —  c 

And,  the  stability  of  growth  following  initiation  is  governed  by  the 
conditions: 


G(a  +  Aa) 

<  V 

stable  growth 

(3-11) 

G(a  +  Aa) 

>GcJ 

unstable  growth 

(3-12) 

G(a  +  Aa) 

~  G  ; 
c 

neutral  growth 

(3-13) 

where  a(>  a^)  is  the  instantaneous  crack  size  for  which  the  growth  stability 
is  to  be  determined. 

3.  Finite  Element  Models 

In  the  foregoing  discussions,  the  available  energy  release  rate  function 
G (a)  was  left  in  general  terms.  For  the  specific  modeling  of  the  transverse 
cracking  and  free  edge  delamination,  G(a)  must  be  appropriately  calculated 
by  the  finite-element  crack-closure  procedures  discussed  earlier.  Details  of 
the  specific  models  are  discussed  in  the  following. 

(a)  Transverse  Cracking  Model  The  finite  element  representation  for  the 
transverse  cracking  is  illustrated  in  Fig.  44.  Under  the  idealized  cracking 
geometry  as  shown  in  Fig.  44(a),  the  growth  path  of  the  crack  is  represented 
by  Fig.  44(b).  The  solution  domain  is  in  the  2-dimensional  x-z  plane. 
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although  the  plane  may  undergo  warping  in  the  y-direction.  Thus,  a  generalized 
plane  strain  formulation  (see,  e.g..  Ref.  [14])  was  adapted  in  the  finite 
element  routine. 

Fig.  45  shows  a  typical  finite  element  mesh  for  the  growth  of  transverse 
cracking.  As  has  been  postulated  previously  that  an  initial  macro-flaw  of 
size  a£  exists  in  the  90°-plies,  and  this  crack  can  be  propagated  to  reach  the 
interface  of  the  face-plies  as  shown  in  the  figure.  Hence,  if  the  laminate  is 
subjected  to  a  uniaxial  tension  represented  by  a  far  field  uniform  displace¬ 
ment  6,  the  energy  release  rate  function  G(a)  may  be  calculated  and  expressed 
by 

G (a)  -  [C  (a)e 2  ]t,  a  >  a  (3-14) 

ex  c 

where  t  is  a  characteristic  length  (e.g.,  thickness  of  a  ply)  and  e^  is  the 
far  field  laminate  strain  induced  by  6.  The  coefficient  function  Cg(a)  is 
calculated  by  imposing  e^  *  1. 

Fig.  46  shows  the  general  behavior  of  Ce(a)  for  the  crack  in  the  lami¬ 
nate  shown.  Note  that  the  crack  size  a  is  limited  by  the  thickness  of  the 
90°-layer,  b.  The  strain  energy  available  increases  initially  with  the  in¬ 
crease  of  cracking  size;  as  the  crack  approaches  the  ply  interface,  the 
available  energy  begins  to  decrease;  the  rate  of  decrease  depends  on  the  rel¬ 
ative  rigidity  of  the  outside  layer.  Generally,  a  maximum  is  located  at  about 
a  ■  0.9b. 

Similarly,  if  the  laminate  is  subjected  to  a  uniform  temperature  drop 
of  At,  and  if  tensile  stress  in  the  90°-layer  is  also  induced,  then  the  calcu¬ 
lated  energy  release  rate  function  for  the  transverse  crack  is  expressed  by 

G(a)  -  [C_(a)AT2]t,  a  >  a  .  (3-15) 

T  C 
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Typical  behavior  of  CT(a)  for  a  transverse  crack  is  similar  to  that  of 
Ce(a),  as  shown  in  Fig.  46.  Again,  C^Ca)  is  generated  by  imposing  AT  ■  1. 

In  practice,  all  laminates  are  subjected  to  thermal  residual  stresses 
before  test;  and  its  effects  on  sub-laminate  cracks  are  not  insignificant. 

The  effect  of  such  residual  stresses  may  be  simply  included  in  the  form  of  a 
uniform  thermal  load,  such  that 

AT  =  T  -  T.  (3-16) 

c  0 

where  Tc  is  the  curing  temperature  and  Tq  the  ambient  testing  temperature. 

Thus,  when  the  laminate  is  loaded  by  unaxial  tension,  a  combined  loading  con¬ 
dition  must  be  considered.  The  calculated  energy  release  rate  function  takes 
the  form 

G(a)  =  [Ce(a)ix2  +  CeT(a)ix(AT)  +  CT(a)(AT)2]t  (3-17) 

where  CeT(a)  is  generated  by  results  for  ex  =  1,  and  AT  =  1. 

Note  that  the  cracking  action  of  the  transverse  cracks  is  primarily  of 
mode-I,  the  opening  mode. 

Accordingly,  for  a  given  laminate  geometry,  say  (+45/0/902)g,  onset  of 
90°-cracking  in  the  90°-layer  can  be  determined  if  the  coefficient  functions 
Ce(a),  CgT(a)  and  C^,(a)  are  calculated.  Then  by  setting  AT  and  to  known 
values  and  by  letting  G(ac)  *  Gc»  Eq.  (3-17)  determines  the  far-field  strain 
ex  at  the  onset  of  cracking.  From  ex>  the  far-field  applied  load  is  then 
obtained;  =  Ex®x- 

(b)  Free  Edge  Delamination  Model  The  finite  element  representation  for  the 
free  edge  delamination  is  illustrated  in  Fig.  47  and  Fig.  48.  For  symmetrically 
stacked  laminates,  edge  delamination  may  occur  along  the  mid-plane  (y-axis), 
or  the  off-mid-plane  interfaces.  The  mid-plane  cracking  is  illustrated  in 
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Fig.  47,  while  the  of f-mid-plane  cracking  is  illustrated  in  Fig.  48.  In  both 
cases,  the  solution  domain  is  the  2-dimensional  y-z  plane,  with  the  applied 
load  in  the  x-dlrection.  Again,  a  generalized  plane  strain  condition  is 
assumed,  if  the  edge  cracking  is  idealized  as  illustrated.  Fig.  49  depicts 
a  typical  finite  element  mesh  for  an  edge  delamination  that  occurs  along  a 
chosen  ply  interface. 

The  edge  crack  begins  with  the  assumed  macroscopic  interface  flaw  of  size 
ac>  and  it  is  propagated  self-similarly  away  from  the  edge  boundary.  If  the 
applied  far  field  load  is  represented  by  the  laminate  strain  e^,  the  calculated 
energy  release  rate  function  again  takes  the  form  of  Eq.  (3-14).  And,  if  a 
combined  thermal  and  mechanical  loading  condition  is  considered,  the  G(a) 
function  takes  the  form  of  Eq.  (3-17). 

A  typical  G(a)  function  for  an  edge  delamination  growth  is  displayed  in 
Fig.  50.  Its  behavior  is  slightly  different  from  that  for  the  transverse 
cracking.  It  is  seen  from  Fig.  50  that  the  available  energy  release  function 
G(a)  increases  initially  with  the  crack  growth  to  reach  a  maximum  value  at 
the  thickness  of  the  layers  containing  the  interface  crack;  beyond  this  point, 
G(a)  remains  at  a  constant  value. 

Hence,  for  a  given  laminate  geometry  and  at  a  chosen  interface  cracking 

site,  the  coefficient  functions  C  (a),  C  m(a)  and  C„(a)  must  first  be  gener- 

e  eT  l 

ated  numerically.  Then  by  given  AT  and  ac  values,  the  applied  far-field  strain 
ex  is  determined  from  (3-17)  by  letting  G(a£)  -*■  G^. 

Note  that  a  mid-plane  cracking  is  primarily  mode-I  action,  while  an  off- 
mid-plane  cracking  is  generally  a  mixed-mode  action.  In  the  latter  case,  G(a) 
is  the  total  energy  release  rate,  consisting  af  all  three  modal  components, 
as  stated  by  Eq.  (3-6). 
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4.  Numerical  Correlation  with  Experimental  Results 

The  AS-3501-06  graphite-epoxy  material  ply  is  assumed  as  a  homogeneous 
medium  having  elastic  orthotropic  property  in  the  principal  material  coordi¬ 
nates.  The  effective  (or  macroscopic)  moduli  of  material  have  been  determined 
in  the  basic  material  characterization  experiments,  see  Table  A.  Specifically, 
these  are  given  as 


140  GPa;  E_  =  E  =  11.1  GPa 
I  z 

GTz  =  GLz  =  5-5  GPa 
VLz  =  °*29 

0. 36xlO-6/°C;  a  =  a  =  28.8xlO-6/°C 
l  z 


(3-18) 


In  addition,  the  average  thickness  of  a  single  ply  is  t  =  0.0136  cm. 

Since  the  laminates  were  fabricated  under  a  curing  temperature  of  177°C, 
and  were  tested  under  the  average  ambient  temperature  of  22°C,  the  equivalent 
thermal  load  in  the  laminates  is  approximately 

AT  =  -  155°C  (3-19) 

The  critical  strain  energy  release  rate  Gc  appropriate  with  the  kinds  of 
cracks  discussed  here  must  be  determined  by  independent  experiments.  Generally 
speaking,  90°-cracking  and  ply  delamination  are  fractures  of  the  matrix,  or 
the  fiber-matrix  interface.  Data  obtained  for  bulk  epoxy  and  epoxy  adhesive 
[30]  under  mode-I  action  range  from  90  to  210  J/m  for  G  ;  and  somewhat  higher 
value  for  G^  under  mixed-mode  actions.  In  a  comprehensive  experimental 

2 

correlation  using  a  T300/934  graphite-epoxy  system,  Law  [31]  found  Gjc~230  J/m  , 

and  a  somewhat  lower  value  for  mixed-mode  G.  .  .  Recent  experiments  by 

t  u*/  c 

Wilkins,  et  al.  [32]  and  Schapery,  et  al.  [33],  have  reported  that  G^c  value 
can  vary  depending  on  the  microscopic  crack  path.  For  example.  Fig.  51 
illustrates  two  situations  where  the  crack  opening  is  in  mode-I,  but  the 
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microscopic  crack-path  in  the  two  cases  is  different.  For  a  graphite-epoxy 

system,  G  value  found  for  cracking  path  parallel  to  the  fiber  direction, 

2 

such  as  shown  in  Fig.  51(a)  ranged  from  80  •*  140  J/m  ,  while  for  cracking  path 

2 

normal  to  the  fibers,  shown  in  Fig.  51(b)  ranged  from  210 r“ 260  J/m  . 

For  the  cases  considered  in  this  report,  the  90°-cracking  action  is 
similar  to  that  illustrated  by  Fig.  51(b);  while,  the  delamination  action  de¬ 
pends  on  the  specific  interface  where  the  crack  occurs.  In  general,  a  mixed¬ 
mode  cracking  occurs  on  off-mid-plane  interface,  and  a  mode-I  cracking  occurs 
on  the  mid-plane  interface  for  symmetric  laminates.  Lacking  the  specific  data 
for  specific  action,  it  is  assumed,  for  purpose  of  illustrating  the  method 
developed  herein,  that  value  for  all  the  different  modal  actions  are  the 
same.  Namely 


G  *  G_  =  G  , 
c  Ic  mixed 


210  J/m2 


(3-20) 


Finally,  the  proper  value  for  the  (macroscopic)  critical  flaw  size,  ac 
must  be  determined.  This  value  is  probably  the  most  ill-defined  quantity, 
and  is  difficult  to  determine  physically.  In  the  experimental  correlation  by 
Law  [31],  ac  for  90°-cracking  was  found  at  about  0.25  mm.  For  delamination, 
the  value  is  slightly  lower.  Thus,  in  the  present  study,  it  is  assumed  that 
ac  is  the  same  for  both  types  of  cracking,  and  its  value  is  given  by 

a  ^  0.25  mm  (3.21) 

c 

With  the  basic  material  values  defined,  the  numerical  correlations  with 
the  experimental  results  are  presented  in  the  next  two  sub-sections. 


*  Mixed-mode  cracking  between  +45°  interface  have  shown  higher  values  than 
this.  But,  carefully  designed  experiments  are  necessary  to  accurately 
determine  these  basic  material  values. 
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(a)  Laminates  under  Static  Tension  Consider  the  experimental  results  for 

the  [+45/0/902 ]g  laminates  tested  under  uniaxial  tension.  Onset  of  transverse 

cracking  occurred  at  131  MPa,  and  onset  of  free-edge  delamination  occurred  at 

262  MPa.  Using  the  basic  material  properties  for  the  material  ply  given  in 

(3-18),  the  finite  element  routine  calculation  for  the  energy  release  rate 

coefficient  functions  are  first  obtained  for  the  90°-cracking  simulation. 

These  coefficient  functions  are  displayed  in  Fig.  52.  It  is  seen  that  the 

three  coefficient  functions  are  all  similar  in  form;  all  have  a  maximum  at 

about  a  =  1.7  t,  which  is  about  85%  of  the  90„-layer. 
m  2. 

It  was  assumed  earlier  that  the  macroscopic  critical  flaw  size  is  about 
a^  0.25  mm,  or  ac'-2t.  But  the  physical  constraint  (the  thickness  of  the 
902  layer)  limits  the  value  of  ac.  Thus,  in  order  to  apply  the  cracking  onset 
criteria  of  Eq.  (3-17),  is  replaced  by  a^.  Hence,  the  equation  governing 
the  onset  of  90°-cracking  is, 

[C  (a )l  2  +  C  _(a  )e  (AT)  +  C_(a  ) (AT)2]t  =  G  (3-22) 

emx  elmx  Tm  c 

Using  the  values  of  C's  from  Fig.  52,  the  value  for  AT  from  (3-19),  the 
value  for  Gc  from  (3-20)  and  t  *  0.136  mm,  Eq.  (3-22)  yields  the  far-field 
laminate  strain  at  onset  of  cracking, 

e  -  3. 15xl0_3. 
x 

Using  E^  from  experiment,  E^  -  43.2  GPa,  the  far-field  laminate  stress 
at  onset  of  cracking  is  given  by 

a  ■  136  MPa 
x 

which  compares  with  the  average  value  of  131  MPa  obtained  from  experiment. 
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As  the  load  Increases,  more  90°-cracks  are  induced  according  to  the 

damage  event  sequence.  Then,  at  some  critical  load,  onset  of  delamination 

occurs.  By  an  edge  stress  analysis,  it  is  inferred  that  the  most  probable 

delamination  sites  are  the  mid-plane  interface  (90°/90°)  and  the  off-mid-plane 

interface  (0°/90°).  Numerical  results  for  the  energy  release  rate  coefficients 

have  been  generated  for  both  cracking  sites.  Fig.  53  displays  the  coefficient 

functions  for  the  mid-plane  delamination  (mode-I  predominant),  while  Fig.  54 

shows  the  coefficient  functions  for  the  off-mid-plane  delamination  (mixed 

mode  I  and  mode  III  predominant) .  By  comparing  these  two  sets  of  curves,  it 

is  determined  that  the  mid-plane  delamination  is  more  energetically  dominant. 

Thus,  using  the  values  of  the  coefficient  functions  at  a  =  a  =  0.25  mm,  and 

c 

AT  *  -  155°C,  G  *  210  J/m 2  in  Eq.  (3-22),  the  onset  of  free  edge  delamination 
c 

along  the  mid-plane  is  obtained  at  the  far-field  laminate  strain, 

i  =  7.25xl0_3. 
x 

Then,  using  the  value  of  of  the  laminate  =  39.6  GPa,  the  far 
field  stress  at  onset  of  delamination  is  predicted  as, 

a  »  289  MPa. 
x 

The  predicted  value  is  about  10%  higher  than  the  experimental  value 
(262  MPa).  The  reason  for  the  discrepancy  is  probably  due  to  the  effect  of 
the  transverse  cracking  on  delamination.  The  interaction  of  the  transverse 
cracking  and  delamination  is  not  accounted  for  in  the  predi^Live  model. 

In  the  experiment,  laminates  [0^/90 2/ ^5 ^-^2  ]g  and  I902^°2^“^52^52^s 
both  showed  transverse  cracking  before  delamlnatlon.  A  numerical  simulation 
of  transverse  cracking  for  both  laminates  has  been  conducted. 
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Fig.  55  shows  the  computed  energy  release  rate  coefficient  functions  for  the 

[0^/90. /45„/-45„]  laminates,  while  Fig.  56  shows  the  same  for  the 
i  i  i  i  s 

[90*/0„/-45«/45» ]  laminates.  It  is  seen  that  the  available  energy  release 
i  i  i  i  s 

rate  for  the  latter  is  much  larger  than  the  former,  because  in  the  former, 

the  90°-layer  is  constrained  on  both  sides.  This  added  constraint  delays  the 

onset  and  the  growth  of  the  transverse  cracking.  Using  the  same  model  as 

before,  the  calculated  onset  stress  for  [0-/90„/45./-45„ J  is 

l  l  Z  is 

a  ~  407  MPa. 
x 

This  value,  however,  is  much  higher  than  the  measured  onset  stress  of 
255  MPa  at  the  first  transverse  crack.  A  closer  examination  of  the  x-radio- 
graph  taken  for  this  series  of  laminates  (3  replicas)  shows  that  significant 
number  of  transverse  cracks  occurred  after  the  load  reached  344  MPa.  On 

the  other  hand,  the  predicted  onset  stress  for  the  [902/02/-452/452 Jg  laminate 
is 

o  =  288  MPa 
x 

which  is  about  20%  higher  than  the  measured  onset  stress  of  234  MPa. 

Finally,  the  predicted  delamination  for  [02/902/452/-452 ]  occurs  at  the 
45/-45  interface,  and  for  [90^ /O2 /— A52 / ^ ^2 ^ s  occurs  at  the  0/90  interface; 
the  predicted  onset  stresses  falls  within  10%  of  the  measured  values. 

(b)  Laminates  under  Static  Compression  Consider  the  laminates  [02/902/452~452 ]g 

and  [90o/0o/-45_/45o ]  for  which  experimental  results  on  delamination  under 
i  i  i  is 

compression  are  available,  see  Table  6.  Onset  stress  is  found  to  be  339.2  MPa, 
for  the  former  and  316  MPa  for  the  latter.  Calculations  using  the  finite 
element  routine  yielded  the  energy  release  rate  coefficient  functions  for 
these  laminates.  In  both  cases,  the  most  energetically  probable 
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delamination  site  is  along  the  +45/-4S  interface,  resulting  in  mixed  mode 
cracking.  The  coefficient  functions  corresponding  to  the  predicted  delamina¬ 
tion  sites  are  displayed  in  Figs.  57  and  58  respectively.  Using  the  values 

2 

of  a  *  a  =  0.25  mm  and  AT  ■  -  155°C,  G  =  210  J/m  and  t  =  2x  (ply  thickness), 
c  c 

the  predicted  far  field  strains  for  at  onset  delamination 

is 

e  =  5 .01xl0~3 

X 

and  for  [90  /0o/45o/-45  ]  is 
i  l  l  is 

e  =  4.65xicf3. 
x 

These  values  translate  to  a ^  =  272  MPa  and  =  240  MPa  respectively. 

It  is  noted  that  the  experimental  values  for  the  former  is  339  MPa,  and  for 

the  latter  316  MPa.  Thus,  the  predicted  onset  of  edge  delamination  under 

compression  is  about  20%  too  low.  The  reason  for  the  rather  poor  correlation 

is  thought  to  be  that  the  compression  specimen  had  a  clear  test  section  length 

of  less  than  2  cm,  so  the  relatively  rigid  end-tabs  exert  a  constraining  force 
normal  to  the  specimen,  which  in  turn  may  suppress  the  onset  and  the  growth  of 
delamination. 

In  the  case  of  the  dispersed  laminates  [0/90/0/90/45/-45/45/-45]g,  a 

45°/-45°  interface  delamination  simulation  would  yield  the  energy  release  rate 

coefficient  functions  identical  to  those  shown  in  Fig.  57,  except  that  t  =  one 

ply  thickness,  instead  of  2  ply  thickness  as  indicated  in  the  figure.  In  this 

case,  the  maximum  available  energy  release  rate  occurs  at  a  ■  t,  which  is  less 

m 

than  the  critical  flaw  size  a^.  Hence,  in  predicting  the  onset  of  delamination, 

the  values  of  the  coefficient  functions  are  taken  at  a  ,  not  a  .  The  pre- 

m  c 

dieted  onset  strain  is  then, 

e  -  8.18xl0~3 
x 
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or,  the  onset  stress, 

a  =  410  MPa. 
x 

The  corresponding  experimental  value  was  380  MPa,  and  the  final  failure 
was  at  425  MPa.  Note  that  at  the  onset  of  delamination  the  specimen  is 
already  near  total  failure.  It  is  not  clear  whether  delamination  or  buck¬ 
ling  action  is  the  cause  of  final  failure. 

In  view  of  the  foregoing  comparisons  between  the  predicted  and  the 
experimental  loads  for  the  various  cracking  events,  it  is  clear  that  the  energy 
method  as  outlined  here  appears  to  be  an  acceptable  analytical  tool  for  de¬ 
scribing  these  cracking  events.  There  are,  clearly,  some  uncertainties  with 
respect  to  the  variable  material  properties.  The  other  limitation  in  the 
method  remains  the  computer  routine  itself.  The  .finite  element  formula¬ 
tion  is  essentially  two-dimensional,  but  most  of  the  cracking  events  observed 
in  the  tests  have  been  three-dimensional  in  nature,  except  a  few  simple 
cases,  such  as  90°-cracks  and  mid-plane  delaraination  in  long  specimens  under 
uniaxial  tension.  This  limitation  precludes  a  prediction  for  those  specimens 
with  implanted  defects.  However,  this  limitation  is  a  numerical  short-coming; 
it  does  not  pertain  to  the  physical  method  of  the  energy  release  rate.  It  is 

conceivable  that  a  three-dimensional,  or  some  quasi-3-dimenslonal  finite  ele¬ 
ment  routine  may  readily  be  developed  to  handle  those  cracks  which  are  three- 
dimensional  in  nature. 
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IV.  FATIGUE  GROWTH  MODEL  BASED  ON  ENERGY 


The  fracture  strain  energy  release  rate  models  developed  in  the  preceding 
section  were  shown  to  be  capable  of  predicting  reasonably  well  the  initiation 
loads  for  multiple  transverse  cracking  and  free  edge  delamination.  Owing  to 
the  2-dimensional  nature  of  the  finite  element  crack  growth  representation,  a 
satisfactory  correlation  between  experiment  and  prediction  for  any  3-dimensional 
delamination  cracking  still  remains  unavailable.  However,  the  energy  model 
does  provide  a  qualitative  means  for  identifying  the  dominant  cracking  modes, 
and  their  dependence  on  the  fracture  energy  release  rate  during  the  course  of 
the  3-dimensional  crack  growth. 

In  this  section,  an  effort  is  made  to  extend  the  concept  of  fracture 
energy  release  rate  to  describe  the  crack  growth  process  for  laminates  under 
fatigue  loads.  Special  emphasis  will  be  placed  on  delamination  growth  under 
constant  amplitude  fatigue.  An  empirical  delamination  growth  model  is  developed 
from  the  experimental  data  that  was  generated  in  early  tests.  The  model  is 
shown  to  apply  well  for  the  tension/ tension  fatigue  case.  For  the  compression/ 
compression  fatigue  case,  the  cracking  is  3-dimensional;  and  the  model  gives 
only  a  qualitative  description. 

A  concept  of  constant  damage  state  is  introduced;  and  by  an  extension  of 
this  concept,  a  tentative  cumulative  damage  model  is  presented.  Details  in 
these  developments  are  described  in  the  following  sub-sections. 

1.  A  Fatigue  Growth  Model  Based  on  Energy  Release  Rate 

For  laminates  undergoing  free  edge  delamination  under  statically  applied 
load,  a  typical  (calculated)  energy  release  rate  curve,  G(a),  is  shown  in 
Fig.  50.  When  the  laminate  is  subjected  to  a  constant  amplitude  fatigue  load. 
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say  the  maximum  available  instantaneous  energy  release  rate  G^(a,o^)  is 
shown  in  Fig.  59.  Assume  that  is  generally  lower  than  the  static  onset 
load  for  delamination;  then  the  laminate  will  not  suffer  any  delamination  dur¬ 
ing  the  first  cycle  of  fatigue  life.  It  is  said  that  the  laminate  possesses 
a  "reserved"  toughness  against  delamination;  and  this  "reserved"  toughness 
will  diminish  with  increased  fatigue  load  a or  increased  fatigue  cycles  N. 

Let  Gc  be  the  material  critical  strain  energy  release  rate  against  delam¬ 
ination,  as  has  been  defined  in  the  static  fracture  case.  Then,  the  quantity 


AG  *  G  —  G£(3,C£) 
err 


(4.1) 


represents  the  "reserved"  energy  release  rate  of  the  laminate  against  delami¬ 
nation  under  the  fatigue  load  a^. 


It  is  further  postulated  that  growth  of  delamination  under  is  described 


by 


a  =  F(AG,N) . 


(4.2) 


Or,  the  increment  of  delamination  Aa  per  increment  of  fatigue  cycle  is 
described  by  the  relation, 


Aa  ,a„n-N 

is  ~  (4G) 


(4.3) 


The  choice  of  (4.3)  is  based  on  the  observation  that,  under  a  small 
fatigue  load  a^,  growth  of  de lamination  is  slow.  From  (4.1)  it  is  seen  that 
for  small  o^,  AG  vG^;  hence  the  left-hand  of  (4.3)  decreases  rapidly  with  N. 
On  the  other  hand,  when  is  large,  or  approaches  the  static  onset  load 
of  delamination  then  AG  <v  0;  and  the  growth  of  delamination  increases 
rapidly  with  N. 

The  exact  form  of  (4.3)  must  be  determined  from  experimental  data. 
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Consider  the  case  of  free  edge  delaminatlon  growth  in  the  [+45/0/902 ]s 
laminates  under  tension  fatigue.  The  averaged  experimental  measurements  of 
delamination  are  plotted  against  log  N  in  Fig.  22,  for  various  fatigue  load 
levels.  Using  the  data  in  this  figure,  a  plot  of  log(Aa/AN)  versus  log(N)  is 
obtained  for  each  of  the  fatigue  loads.  These  plots  are  shown  in  Fig.  60. 
Here 

a  * (a/b) 100%  (4.4) 

where  b  is  the  width  of  the  specimen.  Note  that  a  also  represents  the 
"percentage  delaminatlon." 

It  is  seen  from  Fig.  60  that  for  all  fatigue  load  levels  log(Aa/AN) 
maintains  a  linear  relation  with  log(N).  Then,  the  following  empirical  form 
is  suggested: 

£§■  -  a[k  AG]*0  logN  (4.5) 

AN 

where  AG  is  defined  as  the  "percentage  energy  reserve," 

AG  “  (G  -  GJ/G  .  (4.6) 

c  f  c 

With  (4.5)  and  (4.6),  one  obtains  the  delaminatlon  growth  relation 

a(N)  -  b  (  a[k  AG]"B  logN  dN  (4.7) 

V1 

Strictly  speaking,  AG  is  a  complicated  function;  AG  -  AG(Gc>a,o^)  for 

0  <  a  <  b.  Thus,  (4.7)  is  also  a  complex  integral.  Onset  of  delamination 

is  defined  when  a(N  )  -*  a  . 

c  c 

For  approximation,  one  may  consider  a  »  a and  define 

Gf  -  Gf(am,of)  (4.8) 

where  a  is  the  size  of  a  at  which  Gr  attains  the  maximum  value,  see  Fig.  59. 
m  r 
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Then  AG  becomes  a  constant  and  an  integration  of  (4.7)  from  N  to  N  yields 

m 

a(N)«b  o[N(k  AG)"6  log N  -  y]/ [1  -  31og(k  AG)]  (4.9) 

The  constants  a, 6,  and  k  are  determined  from  experimental  data,  while  y 

is  determined  by  a(N  )  -*•  a  .  It  is  noted  that  the  applied  load  o,  and  the 

mm  I 

structural  geometry  and  their  interaction  with  the  ply  stacking  sequence  are 
implicit  in  AG.  Hence,  the  constants  a,  6  and  k  should  be  primarily  material 
property  dominated.  From  the  data  reduction  plots  in  Fig.  60,  these  three 
constants  are  determined  as 

a  =  0.25;  3  =  0.806;  k  =  12.67. 

and  the  value  of  y  is  y  =  2.35. 

Using  these  values  in  Eq.  (4.9),  and  the  G-curves  generated  for  the 
delamination  cracking  in  the  [+45/0/902]g  laminates  (see  Fig.  53)  growth 
curves  under  various  loads  can  be  calculated  according  to  Eq.  (4.9). 

These  are  shown  along  with  the  experimentally  measured  curves  in  Fig.  61. 

It  is  seen  that  the  growth  model  (4.9)  describes  the  experimental  data  quite 
well. 

Note  that  it  takes  only  fatigue  growth  curves  at  two  load  levels  to 
determine  the  values  for  a,  6  and  k  and  y. 

Next,  consider  the  experimental  data  obtained  from  the  compression  fatigue 
tests  of  the  [02/902/452/-452]g  laminates.  Fig.  29  shows  the  delamination 
growth  curves  for  three  levels  of  fatigue  loads,  at  179  MPa,  221  MPa  and  276  MPa. 
Recall  that  this  series  of  laminates  had  no  edge  delamination  until  a  compressive 
load  of  339  MPa  was  reached.  But  the  predicted  edge  delamination  onset  was 
only  272  MPa.  This  discrepancy  was  attributed  to  the  pinning  effect  of  the 
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end-tabs  of  the  test  specimen,  which  resulted  in  3-dimensional  crack  growth. 

The  finite  element  solution  for  the  available  energy  release  rate  G(a)  however 

was  about  50%  too  high,  because  of  the  2-dimensional  plane  strain  formulation. 

2 

Thus,  if  the  calculated  G,(a  ,or)  is  reduced  by  50%,  and  using  G  -  210  J/m  , 

t  m  t  c 

Eq.  (4.6)  gives  the  values  of  AG  as  0.8,  0.572  and  0.451,  respectively,  for 
of  =*  179  MPa,  221  MPa  and  276  MPa. 

Assume  that  the  delamination  model  under  compression  fatigue  is  also 
given  by  Eq.  (4.5).  Data  reduction  for  this  case  is  displayed  in  Fig.  62.  It 
is  seen  that  the  test  data  generally  follow  a  straight  line  in  the  log(Aa/AN) 
vs.  log(N)  plane.  The  constants  in  Eq.  (4.9)  are  determined  as 

a  -  0.4;  8  *  0.51;  k  ■  35.9;  y  =  4. 

These  values  are  slightly  different  from  those  obtained  for  tension  fatigue; 
under  compression,  the  rate  of  growth  is  faster  than  that  under  tension.  The 
predicted  delamination  growth  by  Eq.  (4.9)  does  not  compare  as  well  with  the 
experiment,  although  it  does  provide  the  general  trend,  see  Fig.  63.  The 
discrepancy  is  believed  due  to  the  inaccuracy  in  the  G(a)  calculation  by  the 
2-D  finite  element. 

2.  The  Concept  of  Constant  Damage  States 

From  the  preceding  discussions  on  delamination  growth  under  both  static 
and  fatigue  loads,  it  appears  that  a  generally  deterministic  relationship  may 
be  found  between  the  applied  load  and  the  state  of  damage.  In  the  event  of 
fatigue  cycling,  this  relationship  varies  smoothly  with  N.  Fig.  64  illustrates 
the  concept  of  constant  damage.  The  vertical  coordinate  of  Fig.  64  represents 
the  static  loading.  When  the  applied  load  reaches  cc,  onset  of  delamination 
occurs.  This  state  of  damage  can  be  independently  produced  by  a  fatigue  load. 
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say  at  a.  <  a  .  Then  under  cj  ,  it  takes  N,  cycles  to  induce  onset  of 
delamination.  Similarly,  the  same  state  of  damage  can  also  be  reached  by  a 
different  fatigue  load,  say  ,  by  Nj  cycles.  Points  defining  the  same  damage 
state  in  the  o^-log(N)  plane  define  a  "constant"  damage  curve.  In  Fig.  64, 
a  curve  representing  20%  delamination  and  a  curve  representing  100%  delami¬ 
nation  are  depicted.  Clearly,  the  100%-curve  corresponds  to  the  final  failure 
of  the  laminate;  and  it  is  the  same  curve  as  the  commonly  generated  fatigue 
S-N  curve. 

The  constant  damage  curves  can  be  constructed  using  experimental  data. 

For  example,  Fig.  65  displays  the  constant  damage  curves  generated  by  tension 
fatigue  for  the  [+45/0/902 ]  laminates.  It  is  seen  that  the  various  damage 
states  (transverse  cracking,  delamination,  etc.)  maintains  a  definite  occur¬ 
rence  sequence  under  both  static  and  fatigue  loadings.  The  shape  of  the 
curves  depicts  the  trend  of  fatigue  degradation. 

Figs.  66  and  67  show  the  constant  damage  curves  in  terms  of  %  delami¬ 
nation  for  [0„/90„/45„/-45„ ]  and  [0/90/0/90/45/-45/45/-45]  -A  laminates, 

L  l  l  l  s  s 

respectively.  These  laminates  were  tested  under  compression  fatigue;  it  is 
seen  that  the  rate  of  fatigue  degradation  is  much  larger  than  that  under 
tension  (see  Fig.  65). 

From  these  constant  damage  curves,  the  growth  of  damage  under  some  spec¬ 
trum  fatigue  loads  may  be  evaluated.  As  has  been  demonstrated  in  the  preced¬ 
ing  sub-section,  a  deterministic  damage  growth  model  based  on  energy  release 
rate  may  be  formulated,  such  as  proposed  by  Eq.  (4.9).  If  such  a  model  correct 
ly  describes  the  damage  growth  processes  in  the  laminate,  it  can  calcu¬ 
late  the  constant  damage  curves  in  the  load-life  plane  witheut  the  necessity 
of  experiment.  The  constant  damage  curves  provide  the  essential  step  in  the 
formulation  of  a  cumulative  damage  model  under  spectrum  fatigue  loading. 
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3.  A  Cumulative  Fatigue  Damage  Model 

Upon  the  construction  of  the  constant  damage  curves  in  the  load  (a^-life  (N) 
plane,  cumulative  damage  in  the  laminate  under  a  given  load  spectrum  may  be  eval¬ 
uated  in  a  manner  similar  to  the  well  known  Miner's  rule.  Consider,  as  an  example, 
the  load  history  shown  in  Fig.  68(a).  When  it  is  applied  to  a  given  laminate,  the 
state  of  damage  at  the  end  of  the  period  of  constant  amplitude  loading  o^ 
is  represented  by  point  A  in  the  load-life  plane.  The  constant  damage  curve 
passing  point  A  is  denoted  as  curve  a^  (say  e.g.,  the  onset  of  de lamina t ion ) . 

At  the  end  of  N^-cycle,  the  fatigue  load  is  raised  to  o2»  Then  the  state 
of  damage  in  the  laminate  at  moves  from  point  A  to  point  B  along  curve  a^. 

Under  and  for  a  period  of  (N2  -  N^)  cycles,  the  damage  state  in  the  laminate 
is  represented  by  point  C;  through  C  passes  the  constant  damage  curve  • 

Clearly,  the  damage  state  represented  by  curve  a2  is  more  extensive  than  the 
damage  state  represented  by  curve  a^.  At  the  end  of  the  N^-cycle,  the  fatigue 
load  is  lowered  to  again;  and  the  corresponding  damage  state  moves  from 
point  C  to  point  D  along  the  constant  damage  curve  a^.  Finally,  at  the  com¬ 
pletion  of  the  loading  history  at  N^,  the  damage  state  moves  from  D  to  E, 
where  the  damage  is  represented  by  curve  a^.  Extension  from  E  toward  the 
failure  envelop  intersects  the  S-N  curve  at  point  F.  EF  gives  the  residual 
life  of  the  laminate  after  it  has  been  subjected  to  the  prescribed  spectrum 
loading  shown  by  Fig.  68(a). 

The  foregoing  cumulative  damage  representation  depends  on  the  ability  to 
accurately  construct  the  constant  damage  curves,  either  experimentally  or 
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Numerically.  In  the  case  of  delamination,  both  aspects  have  been  discussed 
in  this  study.  The  suggested  cumulative  model,  however,  remains  to  be  tested 
against  experiments.  In  view  of  the  understanding  gained  in  this  study,  the 
method  appears  quite  promising. 
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V.  CONCLUDING  REMARKS 

In  this  report,  the  fracture  growth  mechanisms  of  delamination  in  fiber 
reinforced  composite  laminates  have  been  studied.  Experiments  on  various 
types  of  delamination  growth  have  been  conducted  in  order  to  observe  and 
identify  the  dominant  loading,  geometry  and  material  parameters  which  influ¬ 
ence  the  growth  processes  of  delamination.  Predictive  models  based  on  the 
linear  elastic  fracture  mechanics  concept  of  strain  energy  release  rate  are 
then  developed  to  describe  the  various  delamination  processes.  The  energy 
method  has  been  shown  to  apply,  with  varying  degree  of  accuracy,  to  delamination 
cracking  in  laminates  made  of  a  graphite-epoxy  composite  system. 

From  the  numerous  results  obtained  in  the  present  study,  an  assessment  of 
our  understanding  and  ability  to  model  the  complex  process  of  delamination 
is  attempted  here.  The  following  remarks  are  in  order: 

1.  Delamination  stems  from  material  property  mismatch  between  two  layers. 
Under  load,  each  layer  of  the  laminate  would  deform  differently  if 
allowed  to  delaminate.  Bonds  between  layers  prevent  this  from  happening. 
The  bonding  strength  is  derived  primarily  from  the  matrix  phase  which 
is  a  material  property  inherent  to  the  composite  system.  The  driving  force 
for  delamination  derives  from  the  available  strain  energy  release  in  the 
course  of  de lamination  growth.  The  available  energy  near  a  delamination 
crack,  however,  is  a  very  complicated  function  of  loading,  layering 
structural  geometry  and  the  overall  geometry  of  the  laminate.  The  most 
striking  factor  in  the  available  energy  release  is  the  actual  thickness 
of  the  material  layer  adjacent  to  the  delamination  cracking.  This  is  re¬ 
ferred  to  as  the  "thickness  effect"  discussed  in  this  study.  The  thickness 
effect  is  significant  both  in  static  strength  and  in  fatigue  life  of  the 
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laminates.  Accordingly,  the  effect  must  be  considered  In  laminate  design. 
For  all  the  delamination  cases  studied,  except  one,  the  delamination 
growth  front  is  a  2-dimensional  contour  and  the  stress  field  near  the 
contour  cracking  front  is  3-dimensional.  The  exception  is  the  long 
symmetric  laminate  under  tension,  whose  straight  free  edge  undergoes 
delamination  because  of  interlaminar  edge  stresses.  In  this  case,  a 
2-dimensional  generalized  plane  strain  approximation  is  quite  appropriate. 
Consequently,  the  2-dimensional  finite  element  modeling  of  free  edge 
delamination  growth  yielded  exceptionally  good  results  when  compared  to 
experiment.  All  other  cases,  however,  require  3-dimensional  modeling, 
which  so  far  is  not  available.  Nevertheless,  qualitative  analysis  based 
on  the  same  energy  release  method  showed  correct  trend  without  contra¬ 
diction.  Thus,  the  growth  of  delamination  with  2-dimensional  contour 
front  can  be  analyzed  by  the  energy  method,  but  the  quantitative  predic¬ 
tion  must  wait  for  the  successful  development  of  a  3-dimensional,  or  a 
quasi-3-dimensional  numerical  routine. 

The  energy  release  rate  method  provides  a  one-to-one  load-damage  (delami¬ 
nation)  relationship.  It  is  a  deterministic  relationship  based  on  the 
mechanisms  of  crack  growth.  Extension  of  this  model  to  the  growth 
process  under  constant  amplitude  fatigue  loading  is  made  possible  by  the 
concept  of  constant  damage  curves.  Although  fatigue  experiments  con¬ 
ducted  in  this  study  have  been  limited  in  scope,  the  results  obtained 
indicate  the  validity  of  the  constant  damage  concept.  Analytical 

construction  of  the  constant  damage  curves  in  the  load-life  plane  is  also 
possible  if  the  damage  itself  can  be  successfully  modeled  (e.g.,  2-dimen¬ 
sional  vs.  3-dimensional  model  at  the  present  state  of  the  art). 
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The  concept  of  constant  damage  may  be  readily  applicable  in  a 
cumulative  damage  model  for  fatigue  under  spectrum  loading.  All 
this  points  to  a  promising  direction  toward  a  deeper  understanding  of 
failure  in  laminates. 

It  is  noted  that  there  are  several  other  areas  basic  to  the  understand¬ 
ing  of  the  delamination  mechanics,  that  are  left  unexplored.  One  of  which 
is  the  determination  of  the  material  critical  energy  release  rate  G£.  It  is 
conceivable  that  this  property,  when  measured  at  the  macroscopic  scale,  will 
exhibit  a  considerable  scatter.  Compared  to  the  extensive  experiments  on 
the  ultimate  strength  of  composite,  studies  on  the  property  Gc  have  been  rela¬ 
tively  few. 

Another  area  that  needs  more  comprehensive  examination  is  the  fracture 
growth  criterion  under  a  mixed-mode  cracking  action  at  the  macro-scale.  It 
is  apparent  in  the  experiments  that  mode-I  and  mixed-mode  fracture  surfaces 
exhibit  different  characteristics.  Energy  required  to  separate  the  two  dif¬ 
ferent  surfaces  may  be  also  different,  thus  requiring  a  different  set  of  crack 
growth  criteria. 

Finally,  a  concentrated  effort  should  be  made  in  developing  a  3-dimen- 
sional  numerical  routine  which  is  capable  of  modeling  the  delamination  growth, 
similar  to  the  2-dimensional  model  employed  in  this  study.  Such  a  3-dimen- 
sional  routine  is  essential  to  describe  most  other  delamination  cracks  which 
are  3-dimensional  in  nature. 
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Table  1.  Basic  Property  Characterization  Tests 


Type  of  test 


Static  Tension 


Type  of  laminate 


t°l?J 

[9°i2] 

[02/902/452/-452]s 
[902/02/-452/452]s 
[0/90/0/90/ 45/— 45/A5/— 45]c 
[90/0/90/ 0/- 45/45/ -45/45 ]c 
[45/-45/0/90J 

4-  S 

[-45/4 5/90/02]s 


t±A512s 


Static  Compression  [0 


12 


[02/902/452/-4521s 

[902/02/-452/452]s 

[0/90/0/90/45/-45/45/-45] 


[45/-45/0/902]( 

^2s 


#  of  specimens 

2 
3 
2 
T_ 

1 
1 
3 
1 
3 
3 
1 
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Table  2 .  Damage  Growth  Testa  under  Static  Step-Loading 


A  -  Specimens  with  implanted  1/4"  Teflon  strip,  see  Fig.  2. 

B  -  Specimens  t-rith  implanted  1/2"  Teflon  strip,  see  Fig.  2 

C  -  Specimens  with  implanted  1/4"  diameter  Teflon  circle,  see  Fig. 3 
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Table  3.  Damage  Growth  Tests  under  Fatigue  Loading 


Type  of  test 

Type  of  laminate 

ii 

it  of 

specimens 

It  of 

load  levels 

Tension  Fatigue 

[45/-45/0/90J 

Z  s 

9 

5 

[45/-45/0/90J  -  A* 

4-  5 

7 

3 

[45/-45/0/902}g  -  3* 

7 

3 

Compression  Fatigue 

{45/-45/0/902] 

2 

2 

[02/902/452/-452]s 

6 

3 

[02/902/457/"A52Js  -  A* 

5 

3 

1 

[02/902/452/-452]  C* 

4 

2 

[90/0/-45„/'45?] 

5 

3 

[0/90/0/90/ 45/-45/45/  45  ]g 

2 

2 

[0/90/0/90/ 45/-45/45/-45] 

s 

-  A* 

5 

3 

[0/90/0/90/45/— ^5/45/— 45 ] 

s 

-  c* 

4 

i 

2 

A  -  Specimens  with  implanted  1/4"  Teflon  strip 
B  -  Specimens  with  implanted  1/2"  Teflon  strip 
C  -  Specimens  with  implanted  1/4"  diameter  Teflon  circle 
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laminate  lay  up 

avg  ault  ^MPa^ 

E  [GPa] 

X 

mm 

Remarks 

[°12] 

1325.2* 

139.83 

0.291 

*o  without 

max  -  ,  - 

failure 

[9°i2] 

60.0 

11.03 

0.027 

[02/902/452/-«2)s 

524.0 

50.13 

0.283 

[902/02/-452/452]s 

453.0 

51.38 

0.309 

[0/90/0/90/45/-45/45/-45]s 

659.2 

51.99 

0.325 

[90/0/90/0/-45/45/-45/45]s 

673.6 

52.26 

0.304 

[45/-45/0/902ls 

326.1 

39.65 

0.172 

[45/-45/0/902]s* 

430.9 

43.23 

- 

*  new  batch 

[—45/ 45/90/0- ] 

L  S 

930.1 

73.71 

0.317 

[±45]2s 

156.5 

19.10* 

0.69 

*  initial  modulus 

Static  Compression 


laminate  lay  up 

avg  a  lt  [MPa] 

E  [GPa] 

X 

V 

xy 

Remarks 

[°121 

932.2 

136.93 

0.277 

[02/902/452/-452]s 

388.2 

60.33 

0.248 

1 

[902/02/-452/45]s 

346.1 

51.23 

0.292 

[ 0/90/0/ 90/45/-45/45/-45] 

s 

430.9 

53.64 

0.266 

[45/-45/0/90, ] 
l  s 

! 

371.6 

38.13 

0.151 

175.1 


17.30* 


0.60 


*  initial  modulus 
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Table  5 .  Summary  of  Results  from  Damage  Growth  Tests 


Static  Tension 


Remarks 

laminate  lay-up 

E.D. 

[02/9°2/452/-452]s  ! 

255.0 

521.2 

539.0 

+45  cracks  dominant 
from  o  =  427  MPa 

[902/02/-452/452]8 

234.4 

442.6 

541.0 

+45  cracks  from 
a  =  379  MPa 

E.D.  and  +45  cracks 
dominant  from  o  *  448 

[90/0/90/0/-45/45/-45/A5] 

s 

413.6 

- 

650.8 

*  not  clear 

[0/90/ 0/90/ 45/ -45/  45 /— 45  ] 

s 

434.3 

- 

606.7 

*  local 

[45/-45/0/90,] 

131.0 

262.0 

393.0 

E.D.  dominant  from 

o  =  351.6  MPa 

[45/-45/0/90.]  * 
i  s 

131.0 

236.5 

443.3 

*  new  batch 

[45/-45/0/ 90, ]  -  A* 

103.4 

179.2 

435.7 

*  new  batch  implanted 

^  s 

1/4"  teflon  strip 

[45/-45/0/90,]  -  B* 

103.4 

172.4 

434.3 

*  new  batch  implanted 

1/2"  teflon  strip 

[-45/45/90/0,] 

L  S 

820.4* 

- 

935.0 

*  local  T.C.  on  edges 
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Table  6.  Summary  of  Results  from  Damage  Growth  Tests 


Static  Compression 


Laminate  lay-up 

Onset  of  E.D. 
MPa 

a  ,  .  MPa 
ult 

Remarks 

t°l2l 

- 

857 

No  damage  before  failure 

(02/902/452/-452)8 

339.2 

369.5 

Local  +45-cracking 
from  341  MPa 

[02/902/452/-4521s  -  A* 

131.0 

286.0 

Local  +45-cracking 
from  241.3  MPa 

[02/902/452/-452]s  -  C* 

251.6 

330.9 

Delamination  from  edge 
first 

[902/02/-452/452]8 

315.7 

319.2 

[ 0/ 90/0/ 90/45 /-45/45/-45]g 

379.2* 

425.4 

*  Onset  delam.  not  clear 

[0/90/0/90/45/-45/45/-45]  A 

s 

138.0 

382.6 

+45-cracking  from 

310  MPa 

[0/90/0/ 90/ 45/-45/45/ -45]  C 

s 

265.4 

365.4 

+45-cracking  around 
circular  patch 

[+45/0/90,] 

—  1  S 

382.6 

No  delamination  before 
failure 
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Table  7.  Summary  of  Results  from  Damage  Growth  Tests 


Tension-Tension  Fatigue 


laminate  lay-up 


[45/-45/0/90,] 
l  s 

a  _  =  443.3  MPa 
ult 


110.3 


No.  of  No.  of  cycles  T.C.  density  No.  of 

specimens  at  onset  of  at  onset  cycles 

T.C. |  E.D.  of  E.D.  tested 


2,000,000 


1,000,000 


[45/-45/0/90,]  -A 
2  S 

a  =  435.8 
ult 
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Table  8.  Summary  of  Results  from  Damage  Growth  Tests 


Compression-Compression  Fatigue 


laminate  lay-up 

Max  op 
MPa 

No.  of 

specimens 

No.  of  cycles 
at  onset  of 
E.D. 

Max  No. 
of  cycles 

Remarks 

(02/902/452/-452]8 

182.0 

2 

150,000 

1,000,000 

suspended 

a  ,  -  369.5 

ult 

221.3 

2 

3,000 

130,700 

failed 

276.5 

2 

70 

16,975 

ff 

[02/902/452/-452]g-A 

182.7 

1 

- 

257,260 

11 

o  -  *  286.0 

ult 

196.5 

2 

- 

307,000 

suspended 

220.6 

2 

73,430 

failed 

[02/902/452/-452]s-C 

194.4 

2 

10,000 

500,000 

suspended 

a  _  =  330.9 

ult 

222.0 

2 

200 

94,710 

failed 

[0/90/0/90/45/-45/45/-45]g 

276.5 

1 

10,000 

748,320 

f> 

o  , .  =  425.4 
ult 

343.4 

1 

500 

3,100 

II 

[0/90/0/90/45/-45/45/-45]g-A 

219.2 

1 

- 

300,000 

suspended 

a  =  382.6 

ult 

248.9 

2 

- 

143,870 

failed 

276.5 

2 

- 

42,650 

11 

[ 0/90/0/ 90/ 45/-45/45/-45]  -C 

s 

193.7 

2 

5,000 

500,000 

suspended 

a  -  -  365.4 

ult 

219.2 

2 

300 

500,000 

11 

[45/-45/0/90„] 

/  s 

241.3 

1 

- 

100,000 

failed 

a  lfc  =  382.6  MPa 
ult 

304.7 

1 

- 

860 

11 

[902/02/-452/452]s 

180.6 

2 

190,000 

1,000,000 

suspended 

ault  -  319'2 

220.6 

2 

7,000 

77,890 

failed 

274.4 

1 

10,000 

IP 

j 

■I 


19.0 
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Figure  5  Cont' 
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Figure  b  Cont'd:  Load-Sequence  X-radiographs 
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Load-Sequence  X-radiographs  for  [+45/0/902]  Laminate  under  Tension 
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9  Cont'd:  Load-Sequence  X-radiographs  for  [+45/0/90,]  -B  Laminate  under  Tension 
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Figure  9  Cont’d:  Load-Sequence  X- radiographs  for  [+45/0/90,]  -B  Laminate  under  Tension 


Figure  10.  90°-Cracking  and  Edge  Delamination  Growth  Behavior  of 
[+45/0/90-]  laminates  under  Tension 

“*  Z  S 


//-£  6 


U-2  G 


H-2.  6 


II- 2  6 


317.1 


333.3 


357.7 


376.8  [MPa] 


Figure  11.  Load-Sequence  of  X-radlographs  for  [02/902/+452 ]g 

Laminates  under  Step-Compression 
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Figure  12.  Load-Sequence  of  X-radiographs  for  [C>2/902/+452  Js  type  A 
Laminates  under  Step-Compression 
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Figure  14.  Delamination  Growth  vs  Applied  Load  for 
[Q2/902/452/-45^]s  Laminate  Series  under 

Step-Compression 
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Figure  16  Cont’d:  Load-Sequence  X-radiographs  for  [0/90/0/90/45/-45/45/-45] 

type  A  Laminates  under  Step-Compression. 
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Figure  17.  Load-Sequence  X-radiographs  for  [0/90/0/ 90/ 45/— 45/45/— 45 ]s  type  C 

Laminates  under  Step-Compression 
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Figure  18.  Delamination  Growth  vs  Applied  Load  for  the 
[0/90/0/ 90/ 45/— 45/ 45/— 45 1 s  Laminate  Series 
under  Step-Compression 
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Figure  20  Continued 
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Figure  2  3 .  Effect  of  Implanted  Interlaminar  Flaw  on  De  lamination  Growth 
under  Tension  Fatigue  Loading.  [+43/0/ 90.,  ]  . »  type  A  and  tvpe 
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:fect  of  Implanted  Interlaminar  Flaw  on  Delamination  Growth  under  Compression 
Fatigue  Loading.  t0~/90»/45~/-457 ]  ,  type  A  and  type  C. 
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under  276  MPa  Compression  Fatigue  Load 


33  Continued 
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Figure  35.  Life-Sequence  X-radiographs  for  [0/90/Q/90/45/-45/45/-45]g-  C 
Laminate  under  219  MPa  Compression  Fatigue  Load. 


Sure  39.  Effect  of  Ply  Thickness  on  Delamination  Growth  under  Compression  Fatigue  for  type  C  specimens 
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Schematics  of  90°-Cracking 


126 


A 


N ADC-79056-60 


N ADC-79056-60 


Figure  42.  The  Three  Crack  Extension  Modes 


Figure  43,  Finite  Element  Scheme  of  Irwin's  Crack-Closure 
Representation. 
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Direction  of  fibers 


(a)  0°/0°  Delamination  Action 


Figure  51.  Schematics  of  Two  Different  Cracking  Actions  in  Unidirectional  Laminate 


6001 


.Tilt 


id-plane  (mode  I)  De laminat ion  Coefficient  Functions  for  [+45/9/90 
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Figure  58.  Off-Mid-plane  Delamination  (-45/45)  Coefficient  Functions 
for  [902/02/-452/452]8  Laminate  under  Uniaxial  Compression 
and  Uniform  Thermal  Loading. 
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Figure  61.  Comparison  of  the  Edge  Delamination  Growth  Model  with  Experiment  for 
l+45/0/90„]  Laminate  under  Tension  Fatigue. 


Comparison  of  the  Edge  Delamination  Growth  Model  with  Experiment  for 
[0„/90  /45  /-45«]  Laminate  under  Compression  Fatigue. 


65.  Constant  Damage  Curves  Obtained  from  Tests  on  [+45/0/90, 
Laminates  under  Tension  Fatigue. 
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(a)  Schematic  of  A  Spectrum  Loading  History 


Figure  68.  Cumulative  Damage  Based  on  the  Concept  of  Constant  Damage  State. 
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